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Mflss_ and Momentum T urbulent Transport Experiments 
With Confined Swirling Coa xial Jets 


R. Roback 
B..V. Johnson 

SUMMARY 


An experimental study of mixing downstream of swirling coaxial jets 
discharging into an expanded duct was conducted to obtain data for the 
evaluation and improvement of turbulent transport models currently used in 
a variety of computational procedures throughout the combustion community. - 
A combination of laser velocimeter (LV) and laser induced fluorescence 
(LIF) techniques was employed to obtain mean and fluctuating velocity and 
concentration distributions which were used to derive mass and momentum 
turbulent transport parameters currently incorporated into various 
combustor flow models. Flow visualization techniques were also employed to 
determine qualitatively the time dependent characteristics of the flow and 
the scale of turbulence. This study was an extension of the mass and 
momentum turbulent transport experiments with nonswirling. coaxial jets 
which had been conducted in a previous program. 

The flow visualization studies indicated that five major shear regions 
exist in the flow field downstream of the swirling coaxial jets. The first 
four regions were previously observed in the flow field for nonswirling 
coaxial streams. They included: (1) a wake region between the inner and 
annular streams a short distance downstream of the inlet, (2) a large-eddy 
shear region between the inner and annular streams, (3). an annular 
recirculation zone adjacent to the inlet plane, and (4) a reattachment 
region downstream of the annular recirculation zone. A large recirculation 
region along the centerline was also observed for the swirling coaxial 
streams. 


Simultaneous two component velocity measurements were made in the 
radial-axial (r-z) plane, the axial-azimuthal (z-0) plane and the 
radial-azimuthal (r“6) plane with the two color LV system. The velocity 
data pairs were used to determine a momentum transport probability density 
-function (p.d.f.) and average transport rates for each of three measuremen 
planes. A combined LV/L1F system was used to simultaneously measure inner 
stream fluid concentration and one of the three velocity components. Thes. 
concentration/velocity data pairs were used to determine the mass transpori 
p.d.f. and average transport rate in the axial, radial, and azimuthal 
direction. The second central moment (rms fluctuation from the mean), 

skewness, and kurtosis for each mass and momentum transfer data set p.d.f. 
were also determined* 


„ W ? h ® r ® 8ults of the8 e measurements indicated that the largest momentum 
urbulent transport was in the r-z plane. Peak momentum turbulent 
transport rates were approximately the same as those for the. nonswirling 
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flow, condition. The masa turbulent transport process for swirling flow was 
complicated. Mixing occurred in several, steps of axial and radial mass 
transport and was coupled with a large radial mean convective flux. Mixing 
for swirling flow was completed in one-third the length required for 
nonswirling flow. 
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INTRODUCTION 


Background 

Computational procedures -used to describe combustion processes are 
being developed and refined by a number of researchers (e.g., see Ref. 1 
through 5). These computational procedures are used to predict the 
velocity, species concentration, temperature and reaction rate distribution 
within the combustors which in turn are used to determine combustor liner 
heat load, engine performance (combustion efficiency), pollution emissions 
(reaction products) and pattern factor (temperature distribution at turbine 
inlet). Because turbulent flow exists in most combustors of practical 
interest, the calculation procedures usually include mathematical models 
for the turbulent transport of mass (or species), momentum and heat. 
However, the prediction of combustion processes is very sensitive to the 
modeling of the mass and momentum transport processes and improper models 
result in inadequate predictions of combustion efficiency, liner heat load, 
emissions and exit temperature pattern factor. 


The effort in the United States to improve and employ computational 
procedures for the prediction of aircraft gas turbine combustor processes 
is continuing in universities, government laboratories and the aircraft. gas 
turbine industry. As part of the NASA program in this technology field. 
Aero thermal Modeling studies were conducted by three U.S. aircraft gas 
turbine manufacturers (Garrett, General Electric and Pratt & Whitney 
MACA r « ft ? rou ?^ during 1982-1983. These studies were funded through the 
NASA Hot Section Technology (HOST) program and the results are/will be 
available in limited distribution reports (e.g., Ref. 6). The present 
(reported herein) and previous (Ref. 7) experimental studies were conducted 
as part of NASA s aero thermal modeling program* 


The deficiencies in the current computational procedures have been 
attributed to weaknesses in the mathematical models, including the trans- 
port models, and in the numerical methods. One recommendation from a NASA 
workshop on combustion modeling was that the mathematical models used in 
the calculation procedures be validated using experiments specifically 
designed to provide the required input data. (Ref. 1).. A first step in 
this process is the validation of mass and momentum transport models for 
constant density flows.-.. 


The data used to formulate and validate the turbulent transport model s 
have been obtained primarily from velocity and momentum transport measure- 
ments because only "a limited amount of concentration and mass trasnport 
data is available. The mass (species) transport data presently available 
are not sufficient to determine where inadequacies exist in the present 
models or to formulate improvements for the models. One reason for this 
situation is that the methods for. simultaneously obtaining turbulent mass 
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(species) and momentum transport data often have been indirect, requiring 
compromising assumptions. To overcome these limitations, techniques have 
been developed, to measure concentration and velocity simultaneously and 
therefore- obtain mass transport data. This data can be used to evaluate 
and improve combustion oriented turbulent transport -models which include 
scalars such as concentration of species and temperature. 


Several nonintrusive measurement techniques have been -used to obtain 
velocity and concentration simultaneously in recirculating flows (e.g., 

Ref. 8). Raman scattering, marker nephelometry, and laser induced' 
fluorescence (LIF) of a trace material are three examples of techniques 
previously used for obtaining the concentration portion of these measure-, 
ments. In this program, laser induced fluorescence of fluorescein dye as a 
trace element in water was chosen to study the mixing between constant 
density fluids. This technique was chosen for the following- reasons : (1) 
the dye and water are relatively inexpensive, (2) the wavelength required 
to excite the dye is compatible with current LDV equipment, and (3) the 
fluids are convenient to use. This choice restricts the measurement 
technique to the acquisition of constant density transport data. Although 
the combustion process has variable density gases mixing in a reacting 
environment, the mathematical transport models for combustors are expected 
to be based on the turbulent transport phenomena found in constant density 
mixing with modifications for variable density and reacting flows. 

A preliminary effort at UTRC was initiated in 1975 to obtain 
quantitative concentration measurements with fluorescene dye (Paper 28 of 
Ref. 8). .This effort at UTRC was continued in 1978 and made use of 
improved opt ics, data handling capabilities and operative procedures. An 
exploratory study of mixing between confined coaxial jets was reported in 
Ref. 9. The current NASA sponsored study of mass and momentum turbulent 
transport experiments with nonswirling and swirling confined coaxial nets 
was initiated in 1981. The first phase of the study was conducted with 
nonswirling jets and was reported in Ref. 7. This second phase of the 
study, reported herein, was conducted with swirling coaxial jets. The 
current application of the combined LV/LIF measurement techniques along 
with the available data handling procedures provides an opportunity to 
obtain data which can be used to evaluate a number of 'computational methods 
and turbulent transport models. Results from Ref. 7 and the present study 
!? U ® ed evalu *te (1) the widely used two-equation turbulence model, 
U; the Reynolds stress transport model and (3) the probability density 

function formulation for predicting turbulent transport and concentration 
i iuctuations • 


Outline of Present Study 

Turbulent mixing of swirling confined coaxial jets is being studied 
because of its similarity to the combustor situation and thus its value in 
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the mathematical modeling of combustor flow fields. Surveys of previous-, 
experimental studies were presented in Refs. 7 and 10. .The turbulent 
mixing characteristics observed in. the flow field produced by confirm 
swirling coaxial jets are 'applicable to the combustion fluid mixing process 

ch * ra “" l,ticB “• i» sa. turbine combustors end 

furnaces. The coaxial jets provide a method of introducing fuel and air 

,>T bU - t ?- n Cha ® b8r * T he annul * r recirculating flow zones asso- 
ciated with swirling and nonswirlmg coaxial jets expanding in enlarged 

ducts provide the pilot regions usually required to. maintain flame in a 
combustor over a range of operating conditions. 

The Reynolds number (Re = PVd/u) of fluid flowing through various 
sections of aircraft gas turbine combustors vary from 1CT to 1(T and 
!,”! £ ? re %? he fJ° W8 ! rS generally turbulent. The lower Reynolds numbers 
Rpvr.ftiH° r fl K W through ® 00lln 8 holes at engine idle conditions. The higher 
Reynolds numbers occur for How through swirlers or dilution jets at engine 
takeoff conditions. The flow conditions selected for the detailed data 

and^^SoV 11 It 6 ? re8ent 8tud y lead t0 axial Reynolds numbers of 15,900 
and 47,500 for the inner and annular streams, respectively. These Reynolds 

I ZltrZt fac r rs 0f 5 5° 20 than the transitional Reynolds' 

Th» b f S f a ” dare in the ran 8 e occurring in aircraft gas turbines. 

® fore * turbulan t transport phenomena measured in the present experi- 

$a" l ?»rbiue“r *° be typical ° £ tha tr “ 5 P°« Phenomena occurring in 

,f ea ‘“ r “ uhic h have been previously observed (Refs. 7, 10 and 11) 
in the flow fields produced by nonswirling and swirling coaxial jets con- 
fined in an enlarged duct are presented in Fig. 1. For both the^wirling 
and nonswirling inlet conditions, the flow field contains (1) a wake region 
downstream of the annular jet/inner jet interface, (2) a shear layer 
between the inner and annular jets, (3) an annular recirculation cell and 
(4) a reattachment region leading to the fully developed duct flow. The 

C ° at ?i n8 (5) a interline recirculation cell * As Sown 
several iegionsf are relatively complex with interaction between 

The extent of each region pictured in Fig. 1 has been shown to deoend 
upon the dimensions of the coaxial jets and ducts and on the axial and ? 

wuriv utribu ; ions : iheref °"- th * 

iu.V, “r f , i SUall2atl0 " St “‘ iy t0 detennine 'he character of the 
t low and the turbulent transport processes for the experimental 

onf iguration chosen for detailed study. The results of the flow 
visualization study were also used to determine the streamwvse locations 
or obtaining the detailed velocity, concentration and turbulent 
transport rate data. 1 

The major focus of this study was on the acquisition, reduction and 
ly 8 _ 18 _: 0f Velocity * con centra tion, mass turbulent transport rate and 
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momentum turbulent transport rate data at nine axial locations within- the 
duct 'test section. Single component velocity data and inner jet. fluid 
concentration data .were, obtained simultaneously to. determine the .local mass 
(or. scalar) turbulent transport rate. Two velocity components were ob- 
tained simultaneously to determine the. local momentum turbulent transport 
rates. As a result, the concentration and velocity distributions were 
obtained during at least two nonconsecutive data acquisition, runs . The 
data set for each point measurement was analyzed and reduced to obtain the 
mean and three central moments for each probability density function 
(p.d.f.), i.e., the mean values, the rms deviation from the mean, the skew- 
ness nf the p.d.f., and the flatness factor (or kurtosis) of the p.d.f. 

The mean and central moments were also obtained f.or the mass and momentum 
turbulent transport p.d.f.s as well as the velocity and concentration 
p.d.f. s. The reduced results for each data point set are tabulated and 
presented in this report. Graphical presentations of representative data 
are also included to aid in the discussion of the results. The discussion 

,. of these results will be related as -applicable to each region shown in Fig. 

ml* . 
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DESCRIPTION OF TEST APPARATUS AND PROCEDURES 


Test Apparatus 

A schematic Of the test facility used in this experimental program is 
presented in Fig. 2. The principal components of the facility are a water 

®. wate V raa ? fe r and “Bering system, a dye injection system, 
and a test section. For the laser velocimeter tests, the facility was run- 
in a closed, recirculating loop. Water which was at a temperature of 

2 ° C, ^ W ff circulated b y a Pump from the storage tank, through 
metering valves and flow measuring devices to the inner jet and annular 
jet inlets of the inlet plenum. The water in the annular duct and inner 

andV enter ^ d th ! J e8t section, mixed and discharged into the exhaust ducts, 
and ws returned to the storage tank* 

Whenever fluorescein dye was used as a tracer, such as for the flow 
visualization tests and the LV/LIF tests, the facility was operated in a 
single pass mode. The water from the exhaust ducts was discharged into the 
ty sanitary sewer and fresh water replenished the system. For the flow 
visualization tests, dye was added to either the inner jet fluid or the 
annular jet fluid several feet upstream of the entrance to the test section 
to ensure uniform flow of dye into the test section. For the LV/LIF tests 

loif?^d £l °h °l dyS WaS ensured b y. add ing the dye to a mixing chamber 
located a short distance from the dye micrometering valve. The 20 to 40 

dron« r ;n S ?i; e ^ SCr °! 8 the valve was lar8e com P ared to the other pressure 
drops in the system, thus ensuring a constant flow of the dye injected into 

m^nT^f f -r* ? e miXing Chamber was la ^ e to'ensu^ adeqia^e 

mixing and a uniform dye concentration at the entrance of the test section 

‘7*?' i” the dye re.ervoir to keep rte lye 

well stirred and an inline filter was placed in the system to prevent 
clogging of the dye micrometering valve* 

is 8howTin C FiB f 3 he ih 8 V 8 % Cti0n al0ng with the inlet and exhau8t sections 
is shown m Fig. 3. The test section was a 122 mm inside diameter by 1 m 

ong, thm-wall glass tube. When flow visualization and optical exeriments 

lltitllltlf ^ circ " la V Ube8 ’ the water-glass-air interfaces c^n produce 
J 1 Cal . d “ to: c-tion. As shown m Fig. 3, the circular duct test section was 

redice 6 h l* a octangular, glass-walled optical box filled with water to 
reduce beam direction distortion as the laser beams passed from air through 
the glass wall of the duct and into the test section water. A ray tracine 
program was used to determine that the radial displacement of ‘the probe 
volume was less than 0.03 mm and the offset of the measurement direction 
from radial was less than 0.05 deg for radii up to approximately 55 m 

« mn n! at - r t0 - he section entered through an annular duct and a 

mailer inner jet tube. The water then exhausted through the exit duct uo 
over a weir and flowed to the drain. Since the top end^f the e^t dua P 
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containing the weir was open to the -atmosphere, the atmospheric pressure at 
the weir prevented the teat section, from becoming avjrpfessurized. The 
inlet plenum for the annular duct contained three perforated plates to 
produce uniform flow and a honeycomb r cction to remove swirl from the flow. 
The inner jet tube was fed with the same diameter hose for lengths of. over 
300 cm and included a 60-cra length of straight, nonflexible tubing, 
containing a perforated plate, positioned immediately upstream of the inlet 
plenum. 

A detail of the inlet to the test section is shown in Fig. 4. The 
center tube was a 25 mm ID tube with a 7.5 deg half angle flare at the exit 
of the tube. The 59 mm ID annular duct contained a free-vortex swirler 
with 30 deg mean angle flow turning blades whose trailing edges were 
positioned approximately 51 mm from the entrance to the test section. The 
velocity and coordinate system utilized for the measurements is also shown 
in Fig, 4. 

With this test section, measurements were usually made at. a fixed 
axial location over a range of radial locations. The measured parameters 
discussed in this report are presented as functions of the radial position 
r, normalized by the radius of the test section, R , i.e., r/R . For this 
representation, the inner jet tube .extends from r/8 = 0.0 to r/R » 0.20 

and the annular region containing the swirler extends from r/R =°0.25 to 
r/R Q «* 0.48 (see Fig. 4). The inner and annular regions are separated by a 
thin wall which extends from r/R « 0.20 to r/R = 0.25. The region 
between r/R Q s 0.48 and r/R Q = 1°0 is a solid wall which acts like a 
backward facing step in. the°f low field. 


Flow Visualization 

Sketches of the optical arrangements used to obtain flow visualization 
photographs -and motion pictures of the flow pattern in the radial/axial 
plane (r-z) and radial/.azimuthal plane (r-0). are shown in Figs. 5 and 6 
respectively. The light source was an argon ion laser which produced a 1.25 
mm diameter beam and was operated either in single line mode (0.4880 pm 
wavelength) or with all lines operating. The laser beam was passed through 
a glass or plastic round rod which acted as a cylindrical lens and caused 
the beam to diverge in one plane while maintaining a beam thickness of 
approximately 1 mm. The glass rod was positioned perpendicular to the flow 
axis to illuminate the r-z plane along the center line of the test section 
and parallel with the flow axis at selected axial locations to illuminate 
the r-6 plane. Cameras positioned at right angles to the plane being 
illuminated were used to view the flow. In general, the dye concentration 
was increased until the fluorescent light level was just high enough for 
good photographic contrast. When the dye concentration was too high, 
nonuniform light absorption along the light path occurred. 
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LV and_LIF Measurements 


flYgrmy 

The laser veiocimeter (LV) and laser induced fluorescence (LIF) 
measurements were obtained using commercially available components. Some 
electronic components, which were not commercially available when first 
required at UTRC, were designed and fabricated by the UTRC instrumentation 
group. The equipment utilized for each measurement will be.described as 
the technique is discussed. 

The LV measurement system employed in these experiments is sketched in 
Fig. 7 and utilized the two-color LV optics system detailed in Fig. 8. The 
two-color LV concept utilizes the two strong laser lines of an Argon ion 
laser at 0.4880 ym (blue) and 0.5145 urn (green) wavelengths. These two 
colored beams are separated in the optical system and subsequently emitted 
as three beams; a blue beam, a green beam, and a 50-50 mixed blue/green 
(cyan) beam. The three beams are passed through a lens to produce two sets 
of orthogonal interference fringe patterns in one common focal volume each 
having a fringe spacing:. 


d f = X /( 2 sin (0/2)) (1) 


where X is the wavelength of the incident light beam, and 0 is the 
intersection angle between the cyan beam and either the blue or green beam. 
A particle passing through the probe volume will scatter light of both 
colors, blue and green. The light intensity at the photomultiplier is 
modulated by a frequency, f^, corresponding to the particular wavelength 
and fringe spacing and velocity component. This frequency is related to a 
particle velocity component through 


where IK is the velocity component perpendicular to the optical axis and in 
the plane of one set of intersecting colored beams. More detailed 
descriptions of the particular two color laser Doppler velocimetry system 
utilized in this experiment including the frequency shift used to prevent 
flow direction ambiguity may be found in Ref. 12. 

Two Component LV Measurements _ 

For the measurements conducted in this study, the optical system was 
operated in a direct backscattering mode as shown in Fig. 7. The 0.4880 ym 
wavelength beam was used to measure the streamwise or axial velocty 
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component, U. The 0.5145 urn wavelength beam was used to measure (a) the 
azimuthal velocity component, W, when the prole volume was moved 
horizontally across the stream. and (b) the radial velocity component, V, 
when the probe volume was moved vertically. A Bragg cell was used for both 
velocity components to eliminate the flow direction ambiguity. This 
optical subsystem provided signal-to-noise ratios greater than 20 except 
near the test section walls. For the measurements made in this study, the 
nomina? value of the beam intersection angle, 0 , was 9.52 deg. The laser 
beam dia eter was 1.25 mm and the beam Separation at the. 1.94:1. beam 
expander output lens was 53.5 mm. A 310 mm focal length achromatic lens 
was used to focus the beams. With these optical system parameters, the LV 
probe volume was calculated to have dimensions of 0.08 mm diameter, 1.01 mm 
length and contained 28 fringes. 

Besides the sending and receiving optical subsystem, each LV system 
contains other components or subsystems which perform specific functions in 
the flow measurement. These usually include: (1) a scattering particle 

generator or seeder, (2) a traverse system to position the probe volume, 

(3) signal processors, and (4) a data handling subsystem. For the experi- 
ments performed for this study, the particles naturally occurring in the . 
city water supply proved adequate as LV seeds. As indicated in Fig. 7., the 
traverse system consisted of a milling machine base having three directions 
of motion. The range of motion in the streamwise direction was approxi- 
mately 240 mm while the ranges in the vertical and cross stream directions 
were greater than the dimension of the test section. The relative traverse 
position accuracy of this traverse system was approximately 0.1 mm. 

Laser Doppler velocimeter (LDV) signal processors amplify and filter 
the signals from the multiplier tubes, validate the Doppler frequency 
samples, and finally compute the Doppler period which is the reciprocal of 
the. Doppler frequency. - The SCIMETRICS Model 800A signal processors used in 
this study measured the elapsed time for 8 Doppler cycles and recorded the 
pulses from a 125 MHz crystal during the 8 cycle period. The processor 
also measured and recorded pulses for 4 and 5 Doppler cycles, and compared 
them with the 8 cycle result to ensure that the LDV signal. was a valid 
one-particle signal. The integer number transmitted to the data handling 
system is the period of the Doppler frequency in nanoseconds. Two signal 
processors were used in this study (one for each colored light signal). 

Once the LDV signals were processed and accepted, a microcomputer 
data handling system was used to acquire, store and reduce the data. This 
system consisted of (1) a data handling interface (constructed in-house), 
(2) a DEC PDP 10/11 minicomputer with a dual disk operating system, (3) a 
DEC Laboratory Peripheral System (LPS) with an analog to digital (A/D) 
signal converter, and (4) a DEC writer III teletype printer. The LDV Data 
Handling Interface was used to accept only those data points for which the 
two velocity components were obtained within a period of time of 1 msec. 
This time per iod. was considered appropriate for the probe volume length of 
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approximately 1 mm which was used in -this study and for typical velocities 
of 1 m/sec. Data acquisition rate tests conducted under this criteria 
indicated that almost all of the sets of two component velocity data were 
obtained from a single particle moving through the probe volume. A de- 
tailed listing of the equipment employed for the two velocity component LV 
measurements is presented in Table I. 

Combined LV/LIF Measurements 

The tracer dye used for "the LV/LIF measurements was made from fluores- 
cein disodium salt (C£q H 10° jNa„). This dye is used extensively for water 
pollution studies and is availBle from chemical supply houses in powder 
form. Absorption and emission spectra data for fluorescein dye can be 
obtained from Ref. 13. A liquid dye concentrate was produced by dissolving 
2.5 gms of dye powder in 1 liter of water to which was added 1 tablespoon 
of alcohol in order to stabilize the solution. A dilute solution of dye 
made by uniformly diluting 1 ml of concentrate with 3.5 liter of water was 
mixed "inline" with the inner jet fluid in the ratio of 1 part dilute 
solution to 760 parts water. Any variation in dye concentrations at the 
inner jet inlet location can be attributed to this last mixing process. 

The 0.4880 uo wavelength beam of the argon ion laser was used in. the 
LV/LIF experiment both to induce fluorescence of the fluorescein dye for 
the LIF measurement and to scatter light from particles for the LV 
measurements. Fluctuations in the laser beam intensity were monitored in 
bench tests to determine power fluctuations. The peak to peak power drift 
over a 20 minute period was less than 0.5 percent. The signal from the 
photomultiplier was filtered with a 2KHz low pass filter to remove the shot 
noise associated with photomultiplier tubes. The 2KHz filtering was - 
compatible with the typical velocity of 1 m/sec and probe volume length of 
1 mm. The current signal was converted to a voltage, amplified and then 
processed through an A/D voltage converter each time an acceptable LV 
signal was obtained. The LV and LIF data were stored as pairs along with 
the data acquisition time by the Data Handling subsystem. A list of 
additional. equipment used for the LV/LIF measurements is also presented in 
Table I. 


Swirler Design -Evaluation 

Before conducting the flow visualization tests and the detailed mass 
and momentum turbulent transport measurements, a short study was conducted 
to evalute two swirlers designed for use in this program. The eight-bladed 
swirlers had blade shapes which are typical of that used in gas turbines. 
The first swirler was designed for a free vortex flow with a 33 deg swirl 
angle at the mean annulus radius. This swirler produced a flow which 
appeared to separate at the inner tube wall and was deemed unsatisfactory 
for this. program. The second swirler blade shape was designed to produce a 
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30 deg-raean-angle , free-vortex tangential velocity distribution with a 
uniform *xial velocity profile. The second swirler produced a swirling 
annular stream which did not separate from the inner tube wall upstream of 
the inlet plane. In general, the performance of the second swirler was 
satisfactory r.nd the resulting flow patterns were the type desired for this 
study. Therefore, the second swirler was used for the detailed mass and 
momentum turbulent transport measurements and flow visualization teats. A 
detailed description of this swirler is presented .in Fi g , 9 . 
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FLOW VISUALIZATION RESULTS 


The present study was conducted to obtain a data base that can be used 
to evaluate transport models developed for axisymmetric flow computational 
procedures. Since the transport models are based on time-independent 
statistics, it is important that the experimental data used to evaluate 
the models be obtained from statistically steady or stationary flows. 

Also, good axisymmetric flow characteristics are required because data from 
the same radial location are obtained at several azimuthal locations. 
Therefore, flow visualization studies were conducted before data 
acquisition was initiated to determine if the flow was symmetric and 
statistically stationary. The results of the flow visualization studies 
were also used to determine the scale of the turbulent structure of the 
flow within the test section. The structure and scale of the turbulent 
eddies was deduced from the interface between regions of high and low dye 
concentration recorded on high speed motion pictures* 

Flow visualization experiments with swirling flow in the annular 
stream were_ conducted for the following five flow conditions: - 


Flow 

Mean Axial Velocity, m/s 

Flow 

Rate, gpm 

Condition 

Inner Jet, U. 

Annular Jet, U 

a 

Inner Jet 

Annular Jet 

1 

0.52 

1 .66 

6.2 

52.8 

2 

0.27 

1.66 

3.2 

52.8 

3 

2.08 

1.66 

24.6 

52.8 

4 

0.94 

1.51. 

11.1 

48.0 

5 

0.94 

2.87 

11.1 

94.8 


These flow conditions are identicaL to those used in the experiments with 
nonswirling flow in both the inner and annular streams (Ref. 7). 

® ot ion pictures, real-time motion pictures, and fast-action 
U/1000 sec) 35 mm slides were obtained of the flow at selected locations 
in the r-z and r-0 planes. One set of flow visualization tests was 
conducted with a continuous flow of dye added to the inner jet stream. A 
second set of tests was conducted with a pulsed flow of dye added to the 
annular stream to determine (1) the nature of the shear layer between the 
annular stream and the annular recirculation region, and (2) the amount of 
coupling between the large eddy structure in the inner/annular streams and 
in the annular stream/annular recirculation shear layers* 

General features of the flow field observed for all the flow 
conditions except Flow Condition 3 are shown in the sketch in Fig. 1. 

Five distinct flow regions were observed in the flow field. Four of the 
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regions were also observed in the flow field for the nonswirling coaxial 
streams which is also sketched in Fig. 1. They included: (1) a^ake * 

Z'ZlV'cS r UUr acre “"“ * ■*>« distance dolJtream of 

the inlet, (2) a large-eddy shear region between the inner and annular 

streams, (3) an annular recirculation zone adjacent to the inlet plane and 

(4) a reattachment region downstream of the annular recirculation zone’ As 

shown m Fig 1, a large recirculation region along the centerline was’ 

also observed for the swirling coaxial streams. The characteristics of 

Flow Condition i , that condition selected for the detailed data acquisition 

includ!d r i^A in i he J ol i owin 8 Paragraphs. This discussion will also be 
included in Appendix I where the characteristics of the other flow 
conditions will be described. 

Photographs taken for Flow Condition 1 in both the r-z Diane and at 
selected axial locations in the r-0 plane are presented in Fig W The 
Photographs at the top of Fig. 10 .re representative of the flow ' 
characteristics as seen in the r-z plane when dye is injected both in the 
inner stream and in the annular stream. The photograph taken when dve is 
injected into the inner stream (top left) show's tt£ a hi‘gh inner jtl 
stream concentration persists for an axial distance of approximated 50 mm 

riehn P J oto f rap J taken when dye is injected into the annular streak (top 
right), dye has been entrained near the centerline between an axial P 
distance of approximately 60 mm to 170 mm indicating the presence of a re- 
circulation cell near the centerline. The existence of this recirculation 

a^t~Xr r %o?rTJ “ *J* .ppear.no sin, 

Imost 150 urn. Both photographs show large-scale eddy structure in the 
shear layer between the jets starting at axial distances of approxLatelv 
eSdU.n eXte ^ ln8 downstream to at least 175 mm. The scale of Sse 7 

clearly def ined.^he “ e • 

Z'lVl™ f low^condi t ion . 1 ^ obs-^ 1 ^^ ^rLliLung^ 

radial h scSle t (i r ? P ™) t of^h° f hh* ^ the r_6 plane at 2 “ 25 “m shows the. 
racial scale 1= 6 mm) of the eddy structure at the upstream end of tba 

shear region. At z - 51 mm, the scale of the large eddy “ruc^urf El 

increased significantly (^15 to 20 mm). The scale does not all h ! 

increase appreciably fiom z - 100 mm^o z - 200 ^ It shllnll tl *°v 

graphs in both the r-a and r-0 planes, the l.rgeTddy .trSnre'is Lf ' 

axisymmetnc or azimuthally periodic for z > 50 mm. 7 ructure 18 not 

8ecti^ 0 ^f? Ph -tw ta i en f t . 8elected axial locations downstream of the test 

These nW *1 l ? jected in the annular stream are shown in Fig? 11 

These photographs show the structure of the eddies in the r-O^lane. It 
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z ■ 13 nun, the flow on both the ID and OD of the annular stream show small 
scale structure with the type of wall eddies expected in an annular duct. 
At z ■ 25 mm, the scale of the eddies in the annular jet/inner jet inter- 
face increasd to 6 to 20 mm. Within the recirculation cell (38 mm < z < 
180 mm) the scale of the eddies varied from 5 to 25 mm or up to 20 percent 
of the test section duct diameter. Downstream of the recirculation cell 
(z > 180 mm), the flow initially had large lobes (e.g., z ■ 203 mm) which 
evolved into a vortex swirl pattern (z = 406 mm). There was little radial 
mixing downstream of the recirculation cell. 

The high speed motion pictures also provided insight into the mixing 
process and showed two major mixing regions: (1) at the interface between 
the inner jet and the centerline recirculation zone (z * 30 to 60 mm) and 
(2) in the interface between the inner jet/annular jet streams (30 mm < z 
<90 mm and 0.3 < r/R^ < 1.0). The mixing at the inner stream/recircula- 
tion cell interface diluted the inner jet concentration and the resulting 
mixture of fluid from the recirculation cell and the inner jet was 
entrained by the swirling annular stream. 
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FOREWORD TO PRESENTATION OF RESULTS 

r? uct ion „r 

data obtained in this study in addition^* t0 b ! detennined from the 

velocity components and'^entm^Tus^alt SSiSj , „ 

(1) parameters which can be used to characterise th k included 

functions of the velocity components ehl J the probability density 

port rates and the momentum transport Jates and”(2f Jh» 8 ' maSa tra ”*' 

coefficients for the transport processes! Varl °“ a cortelatio ” 

able 

52 tatS'rlST'S Zt T” titie5 hava baa " oalculated'and are 

and three “ent^l ZeJis o£ the preae,,ted include the «* » value 

density functions ?!! the velocity and concentration probability 

or flatness factor)/the mean^alue^nd^h 111 ^ mean » skewness and kurtosis 
and momentum turbulent transport rate nmh^M ^ nt J al moments of the mass 
the correlation_coefficienta P for li P rob ability density functions, and 
rates. &tl0n ~ coe£ flcienta for the mass and momentum -turbulent transport 

for n« SSi^iXSta.'S ?L e 8 n ted in ‘? ia rap °« ««• 

=£? = =JS3 2 rs 

location or change in measured parsmeter was made- i e Zi* aXlal 
or concentration. The number of velocitv/vei^i!’ 1 ,e ** velocity component 
tion data pairs which were acnniraH a • ^°o it y or velocity/concentra- 
was either 250 , 500 or 1000 dene-id' Uring ea £ h sin 8 le point measurement 
traversing the V°°° de P endin 8 “Pon the number of particles 

on floppy disks! Thirdata*wa! UrX h 8 ^ ? c< l ui8ition > all ^ata was stored - 

lated parameters llj.il Jn the ,£Z£ n ‘ lr r,d “* d ' 10 obtai “ ‘*> a «lcu- 
pairs actually used in Him a *• ^ j ? paragraph . The number of data 

500 or 1000 data pairs acquired bec UCtl °5 T” U ? Ually fess than the 250 , 
data reduction whenever spurious dat * 186 * ta pairs were eliminated during 
defined a. data n“c"n”i ~ ” ^"Tt^t rKrf Z‘ ~ 

pas? SsSf w 

are presented'iil^hu'Je^rirub^urfoJm 10 ^' 18 ^ J**f ac 1 uii| rtion run 
presented in graphical form to aid in the di^cu^Hnlf rte «^Us"V l8 ° 
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listing of the run numbers from which data was utilized for the tables and 
figures presented in this report is presented in Table II. A listing of - 
numbers assigned to the figures on which results are displayed in this 
report is presented in Table III. The tabulated parameters for each data 
point set are tabulated in Tables IV-XX where the term XX denotes the run 
number . 

The results are presented and discussed in the following order. The 
mean and fluctuating velocity and concentration results are presented 
first, the turbulent mass and momentum transport rates and correlations 
second, the higher moments of velocity and concentrations, and the higher 
moments of the turbulent transport. rates .are presented, in later sections. 
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DISCUSSION OF MEAN AND FLUCTUATING VELOCITY AND CONCENTRATION RESULTS 


Velocity Results 

Mean and fluctuating velocity profiles were obtained at nine axial 
stations as part of both the mass and momentum turbulent transport measure- 
ments. Consequently, each velocity profile is comprised of data obtained 
in two to four radial surveys (two or more runs) through the center of the 
test section. The coordinate system employed for this study is presented 
in Fig. 4. The results will be discussed in relation to -the shear regions 
shown in Fig. 1 and will be compared with the results of the study made for 
nonswirling coaxial jets presented in Ref. 7. 

Mean Axial Velocity 

The mean axial velocity profile along the centerline of the test - 
section is shown in Fig. 12 and is compared with the profile obtained for 
the nonswirling flow. As can be seen, the velocity profile for the 
swirling flow is .significantly different from that observed with the 
nonswirling flow. For the nonswirling flow, the centerline velocity was 
(1) initially decreased as the inner jet fluid momentum is transported into 
the wake between the inner and annular streams, (2) increased as the inner 
jet fluid is accelerated by the annular jet until z * 200 tom, and (3) 
finally decreased as the duct velocity profile begins to approach that for 
fully developed duct flow. For the swirling flow, the velocity along the 
centerline (1) decreased rapidly, (2) passed through zero at z = 40 mm, 

(3) remained negative until approximately z * 167 mm and (4) gradually 
increased to the average value for fully developed duct flow. The two 
stagnation points at z « 40 mm and -z = 167 mm mark the extent of a recircu- 
lation zone near the centerline of the test section. The length of this 
zone, J.27 mm, agrees with the length observed in the flow visualization 
tests for this flow condition. 

The mean axial velocity profiles are presented in Fig. 13. The 
velocity profile at the measurement location closest to the inlet, z = 5 mm 
was similar to that observed with the nonswirling flow for z » 13 mm. For 
both cases, the peak velocity in the inner jet was approximately one-half 
the peak velocity in the annular jet and there was a wake region in the 
vicinity of r/R = 0.25, the radial location of the thin wall which 
separates the inner and annular jets. 

For the horizontal traverses, the shaded symbols shown in the figure 
were obtained on the side of the centerline adjacent to the LV optics while 
the open symbols were obtained on the far side. For the vertical 
traverses, the open symbols were obtained from the upper half of the 
traverse and the shaded symbols were obtained from the lower half. The 
velocity profiles are axisymmetric since data taken on either side of the 
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centerline generally fell on top of one another. However, there is some 
discrepancy, in the annular region for z - 5 mm. This discrepancy may be 
attributed to the swirl vane wakes and the secondary flow which forms 
adjacent to the swirl vanes and persists in the flow for a short distance 
downstream from the test section entrance. This behavior will be discussed 
more fully later (Appendix II) when the results are presented of azimuthal 
and radial traverses of the flow field a short distance (z » 5 mm) 
downstream of the trailing edge of the swirler. 

The changes in the axial velocity profiles from z = 5 mm to succesive 
downstream locations document the development of the various shear and re- 
circulation regions within the test section. The velocity of the inner jet 
decreased with axial distance and became negative indicating a reversal of 
flow near the centerline region. Note at z = 25 mm, the -velocities were 
negative in the wake region at 0.18 < r/R < 0.31 indicating an "s" shaped 
path for the fluid in this region. The f?ow visualization photographs and. 
motion pictures showed significant mixing in this region. A second flow 
reversal in this region occurred between z ■ -152 mm and z = 203 mm after 
which the velocity gradually increased and flattened out throughout the 
rest of the duct length. 

Negative axial velocities near the peripheral wall were also observed 
for axial locations up through z =51 mm indicating back flow or a 
recirculation region near the outside wall. Note that the peak negative 
velocity near the wall was approximately 0.65 m/s and was measured at an 
axial location of z = 51 mm. This negative velocity is approximately 47 
percent of the peak streamwise axial velocity at this location and is 
almost three times the percentage obtained for the nonswirling case. Also, 
the downstream end of this outside wall recirculation zone is located 
between z * 51 and z = 102 mm whereas the downstream end of the 
recirculation region for the nonswirling case was located at z = 254 mm. 
Therefore, the recirculating zone near the outside wall is shorter for the 
swirling case and the recirculating flow within the zone is stronger. One 
other difference which 3hould be noted is that once reattachment of the 
annular jet occurred the axial velocity profile tended to flatten with 
momentum transport from the outside inward for the swirling flow whereas it 
tended to flatten with momentum transport from the inside outward for the 
nonswirling flow. 

Mean Radial Velocity 

The profiles of the mean radial velocity are presented in Fig. 14. 
These profiles were generally axisymmetric except in the annular region at 
an axial loction of 5 mm. This behavior has been attributed to swirl vane 
wakes and/or secondary flows which form adjacent to the swirl vanes and 
persist in the flow for a short distance downstream of the test section 
inlet (see Appendix II). The flow showed significant radially oucward 
movement (positive velocities) for z < 51 mm. At z = 51 mm, the -dial 
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Si 0 ®}? ia fl PP^i”«tely 75 percent of the axial velocity indicating that 
low was moving outward at an angle of approximately 35 degrees. 8 ’ At 
* " ” ”»• th «e radially inward movement of the flow -near the wall 
r/R o l 0,7 ® 8 ® in “Seating the presence of the recirculation zone which 
was observed in the flow visualization tests and which was also indicated 

102 ^ I ri 0Clti r at this location * Between z = 51 and z * 

ram, the flow attached to the wall and was directed radially inward as 

the flow .a. deflected by the wall. The radial velocity profile Hastened 
out and approached aero aa the duct flow became fully developed (a >305 

Mean Azimuthal Velocity 

The mean azimuthal velocity profiles are presented in Fig. 15. The 

onp m »n hal i Pr0flle8 I®” 8 ! nerally ^isy^etric and well behaved. However, 
T y ? CCUrred l *., the velocit y Profile near the inlet (z = 5 mm) 
where the velocity profile on one side of the centerline 0.4 < r/R <05 

ZirT ?*? ° n ?! ° ther 8ide * Similar behavior waa a l a ° notgd in'the 
axial and radial profiles for this location (see Figs. 13 and 14) and has 

d \riiT rl vane wakes and/or 8econdary fiows in the 

test ^t SWirUn8 fl ° W azimuthal velocity profile near the entrance to the 
TWe Wa - f 8 . ex P® ct ® d fr °m the swirler design characteristics. 

There w^s no swirl in the inner. jet. region and the annular profile was 
essentially that which would be predicted from the swirler design. The 
azimuthal velocity decreased over the first 51 mm after which an 

“! Uy Unif0ra profile was maintained throughout the remaining length 


Fl uctuating Velocities 

Fluctuating axial, radial, and azimuthal velocity profiles are ore- 
sented m Figs. 16, 17 and 18 respectively. The profiles for the swirling 

t0 th ° Se observed for the nonswirling jets in that they 
could be related to the developing shear layers. Near the entrance to the 
est section (z = 5 mm), the peak fluctuating velocities occurred in the 
wake region separating the inner and annular flows (r/R « 0.25) and in the 
°“f side the annular jet (r/R, - 0.55). The !uteneitl4 of the 
fwi T 8 v ®}? c * tlas inc reased and decreased with the development of the 

ainila^ierand'fi 6 ^ 66 ” th ? a " d annular i ets * < 2 > between the 

annular jet and the recirculation zone near the outside wall and m in 

reatt.chment region. The magnitude of the fluctuation, “re generally 

!wir iL now th‘r iU " 8 - l0W tha ;, th,t f ° r the non ®vtrling flow" For the 
mTwhi? 8 J, h m0re VI 8° r0U8 fluctuations occurred within the first 51 
mm while the fluctuations increased to z « 152 mm for the nonswirling 

between sTram an^inJ 108 tbe fl J* ctuations dampened relatively rapidly 

etween 51 mm and 102 mm which might be attributed to the stabilizing ? * 
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effect _ of the uniform azimuthal. velocity profile. The fluctuating radial 
? nd aaim Y t * ial velocity profiles showed trends similar to the axial profiles 
but the intensity of the axial velocity fluctuation was generally greater. 
As the local shear rate decreased (z > 152 mm), the ratios of fluctuation. 
v J n a " d te ? ded toward 1.0. However, for those axial locations- 

where the fluctuations were more intense, the axial fluctuations were 1.5 
o 2 times greater than either the radial or azimuthal fluctuations 
indicating nonisotropic turbulence. 

Concentration Results 

For the concentration measurements, a small amount of fluorescein dye 
was added to the inner jet fluid. As the dye in the inner stream mixed 
with the annular stream, the local concentration of inner jet fluid was 
determined by measuring the intensity of light emitted by the fluorescing 
k u 6 *Tir 6 "Shb intensity is proportional to the concentration of dye in 
l \ l TObe Y olume \ The local concentration in the fluid, f, was defined 
to be the ratio of light emitted locally to the light emitted at the 
centeriine locations at an axial location of. 13 mm where f * 1.0 by 
definition. The location z = 13 mm and r/R = 0 was chosen as a reference 
ocation because: (1) optical interference with the upstream walls occurred 

£ . C ° lle ^ tl0n system at locations closer to test section inlet 
and (2) the inner jet fluid had not begun to mix at this location. During 
a particular run, the light intensity at the reference location was 
measured at the start of the run, after each set of approximately 5 data 
points had been taken, and. again at the end of the run. In this way, any 
variation m the light emitted at the baseline location could be accounted 
darin 8 the course of the run. In the discussion of the experimental 
results which follows, the symbol "f" and/or the term "concentration” refer 
to the ratio of the light emitted locally to the light emitted at the 
baseline location. 


Mean Co ncentration 

The axial variation of the concentration along the test section 
centerline is presented in Fig. 19 along with the variation obtained for 
the non-swirlmg flow (Ref. 7). Again, the results presented in this 
figure point out the significant differences between the Swirling and 
onswirlmg flows. For the swirling flow, the mean centerline 
-°n e S at T decrea f ed ra Pi d ly from the reference measurement location (z 

0 lllTll Tt fnS r ° a ' ^ the ® a8S averaged flow concentration level of f = 
0.104 at z 100 mm. For the nonswirling flow, the concentration alone the 
centerline decreased slowly to 75 mm and then fairly rapidly from z *• 100 
mm to z = 200 mm and finally approached the mass averaged flow 
concentration level at z = 356 mm. Mixing along the centerline Of the test 
section for swirling flow was therefore completed in one-third the length 
required for nonswirling flows. 8 ” 
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preflenSd in Pi^ IS Ji Proftlea at yarioua axial stations are 
section (a*- n 8 l,? 0 ^t, At an axial location near the entrance to the teat 
i*. 13 “"?!* tha concentration profile shows the expected large 

concentration gradient between the inner jet and the annular jet. Note 
that at z ■* 25 ram and 0.18 < r/R < 0.31, the location where 2 wciSSa- 
tion .zone is indicated by the mean axial velocity profile (Fig. 13) the 
concentration profile varies from f - 0.8 to 0.3. The clnce^ration’ 
profile rapidly flattened out with axial distance until z = 1Q2 mm where 
the concentration across the duct became essentially equal to the average 
duct concentration of f - 0.104 for the remaining length of test section 

reLI e8 hi t8 8 ^ 0Wn in ? i8 ‘ 20 indicate that th e concentration profiles were 
reasonably axisymmetnc and repeatable. P we e 

Fluctuating Concentrationa 

a xi fll T wT Centrati °u fluctuationa » are presented in Fig. 21. At an 

SnnrLJ !- 10n r? ear th ? entrance to the test section (z = 13 mm), the peak 
concentration fluctuation occurred in the interface region between the 

= 25 r n^ nd th nnUla J JetS where . a lar ® e concentration gradient exists. At z 
25 mm, the peak concentration fluctuation still occurred in this 
interface region (0.15 < r/R < 0.35) but increased ability alL o.c, 
along the centerline, r/R =0.0. The fluctuations in the centerline 
region continued to increase between z = 25 mm and z =51 mm with peak 
activity shifting toward the centerline. There was a rap^d"^^^ 
activity between z =51 mm and z = 1Q2 mm as the fluctuations dampened out 
from inside to outside. Between z = 102 mm and z = 203 

wall rIIion re ThHl^V lat <f with ? nl y sli 8 ht activity near the outside 
11 region. The fluctuation for swirling flow were essentially dampened 

flow 7 a PP r0Xlmatel y one-third the length required for the nonsSirli!* 


occurred 
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DISCUSSION OF TURBULENT TRANSPORT RESULTS 


Momentum Transport 

t* 16 momentum turbulent transport rates and correlation 

nl 7 ei £ T thlie T-' « pSSiSu 

22 t ' ,I0 “8 h 27 • Ihe discussion of these profiles will be related to 
iX nT^h P [r“ ted 10 FiS - 4 “" d “• presented 

Radial-Axial Plane (r-gl 

^ he “omentum turbulent transport rate -distributions obtained in the 
rllfltiH * 5 ! pr ® 8e ? ted ln 22 • This turbulent transport, uv, is 

transport rates were either positive (outward) or negative (inward) at 
various locations in the flow and were consistent with the local velocitv 
gradients shown in Figs. 13 and 14. Near the entrance of the Jest section 
2 T the oomentum turbulent transport was positive for r/R <02 (th« 

Sri ?h ?i ° n) “! P °' itlve UP ~ 

location of z * 25 mm. This momentum flux can be attributed to the shear 
of the inner jet fluid on the wake region downstream of the interface 
between the inner and annular jets. At z = 5 aim, and for 0.2 < r/R <04 
the momentum flux was negative because momentum was being transported f^om 
the annular jet fluid into the wake region. Near r/R - 0 . 5 ? tJe momentum 

be«.eh 8 thT^a a r end )h “ tr ? bu “‘ i “ «>• layer which develop, 

eeween tne annular jet and the recirculation zone near the outside wall 

had approximately 'doubled? = 25 momentum turbulent transport rat^s 


• w,- T *“\ fl ? W visualization tests and the axial velocity measurements 
indited that a recirculation none existed near the teat !ecUon center- 
line for axial location. 40 am, < a < 167 aim. The result, shorn! in Fia ’2 
indicate that there is little momentum transport near the centerline, r/R 
< 0.4, for axial locations, z = 51 mm and z = 102 mm, the upstream end o? 

for r/R e < 1 0 Il 4 wh^h - The f e Was a weak ne 8 at ive momentum flux 

, ~. 152 111111 whlch 18 cl08e t0 the downstream edge of the 

the o^tsiL r wan CU r/R 10 > ITl Th ®. “^tum turbulent transport rate near 
uue uutsiae wail! r/R > O.o was still positive at z » si mm k..#. 

i 

z e < 8 203 mGnt8 t^ e '^ eW ^° wt ^ 8 ^ a ^ Za ^i°^ "ests^^lt^xiafioc^tiins^M mm < 
£ith 2heT\^* m ° men f tum turbulent transport rate was negative (consistent 
with the flAial velocity gradient) and was essentially zero for z > 305 mm 
In comparison, the momentum transport in the r-z plane for nonswirliSf 
flow was generally positive and_signif icant transport was still observed at 
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an axial location z - 305 mm* The peak rates for both swirling and 
nonswirling flow were of the same order of magnitude ( 0.6 m Z / 8 2 ) but the 

peak for the swirling flow occurred at an axial location one-fifth to 

one-fourth of that observed for the nonswirling flow. 

tr«nn!nrt c ® r re!La t ion coefficients, R , obtained for the momentum turbulent 
transport data are presented in Fig. u 53. As can be seen from the figure, 
the correlation coefficients generally followed the same trend as the 
momentum turbulent transport rate profiles. Peak correlation coefficients 
occurred at essentially the same locations where the peak transport rates .. 
occurred. The absolute values of the peak momentum transport correlation 
coefficients were between 0.35 and 0.50 which are in the range previously 
measured for turbulent free shear momentum transport (Ref. 14). Note that 
the momentum turbulent transport coefficients for z = 51 mm and 0.4 < r/R 

< J’ 6 fT® “°* 25 t0 ~ 0 * 30 rather than the values of 0.4 to 0.5 measured in 0 
other high shear regions. Note also that for z > 305 mm the correlation 
coefficients were low and scattered ahout zero. 

Axial-Azimuthal Plane (z-fi) 

The axial-azimuthal momentum turbulent transport rates, obtained are 

(l/I^au/a/^u/s’ T 5 iS tu y bulent transport is related to the strain rate 
U/r) au/30 + 3W/3z. For axisymmetnc flow, 3U/30-, is zero. The trans- 
port rates in this measurement plane were consistent with the trends shown 
with the corresponding velocity profiles but the peak rates were about half 
of those observed in the r-z plane. The momentum turbulent transport rate 
was generally Positive for z < 102 mm especially in the annular jet region, 

0,25 y .0 < °* 60, For 102 mm < z < 203 mm, the momentum transport rate 
was positive near the centerline region, r/R < 0.4 and negative near the 
outside wall, r/R Q > 0.6. The momentum turbulent transport rate in the z-0 
measurement plane became essentially zero for z > 305 mm and r/R >02 

Iectio^ 8 r/R We < k 0?2! atiVe transp0rt close t0 the centerline of tKe test* 
o . 

The momentum transport correlation coefficient, R , obtained in this 
measurement plane is shown in Fig. 25. Again, the collation coefficient 
profiles generally followed the same trend as the momentum transport rate 

‘ "Jr"' the pC * k values those obtained for“he 

r-z Plane. There appears to be more scatter in the R correlation coeffi- 
cients than those for R . The flow visualization motion pictures showed 
more low frequency "wanting" of the flow inside the downstream of the 

Cf*Tl tffl 1 mil of i ^ i. _ _ . i * 


UW 


central recirculation zone that may cause this scatter in R 
Radial -Azimuthal Plane (r-6) 

The momentum turbulent transport rates in the radial-azimuthal plane 
were obtained by aligning the laser beams in the test section parallel with 
the axis of symmetry. For these measurements, the LV laser beams were 
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passed through s 1200 nun focal length lens and entered the test section 
through a window at the downstream end of the teat section. Light from the 
scattering particles in the probe volume was collected by a set of 
receiving optics which was aligned perpendicular to the test section axis 
of symmetry^ With this arrangement, data was obtained at three axial 
locations: z » 25 mm (Runs 47, 48 and 49), z - 50 ram (Runs 46, 50, 51, 52 

and 53) and z = 102 ram (Runs 54, 55, 56). Runs 46, 47, 50, 51 and 54 were 
obtained with the LV optical axis parallel to the test section axis of 
symmetry. Because of some optical axis constraints with this particular 
optical orientation, data could only be obtained in the region r/R < 0.5. 
In order to obtain data in the region r/R >0.5, the optical axis°was 
tilted upward or downward by approximately one degree. This tilted optical 
orientation allowed data to be obtained closer to the outer wall, Runs 48, 

^2, 53, 55 and 56. This data should be corrected using the result of 
the momentum turbulent transport data obtained for the r-z plane. However, 
since this correction is small (estimated to be within 2 percent of the 
peak turbulent transport rate measured) and since the details (i.e., higher 
moments of the vw probability density function would not be available with 
the correction, uncorrected data is presented herein. It is believed that 
the uncorrected data provide an adequate representation of the momentum 
turbulent transport process in the r-9 plane. 

The momentum turbulent transport rates obtained for the r-0 
measurement plane are presented in. Fig. 26. This turbulent transport is 
usually related to the strain rate (1/r) 3 V./ 30 + 3W/3r - W/r. For 
axisymmetric flow, the azimuthal variation of V is zero. This data was 
obtained for three axial locations, z - 25, 51 and 102 mm. At z = 25 mm, 
there was a negative transport rate in the interface region between the 
inner and annular jets and a positive transport rate in the shear region 
between. the annular jet and the recirculation zone near the outside wall. 
The results obtained for the r-0 plane were similar to those obtained for 
the r-z and z-0 planes in that peak momentum transport occurred at an axial 
location z * 25 mm and the rate tended to zero with increasing axial dis- 
tance. The peak rate observed for the r-0 plane was approximately equal to 
the peak obtained for the z-e plane Ov 0.3 m /s ). 

The profiles of the correlation coefficient obtained in this plane, 
shown in Fig. 27 and generally followed the same trends as the 
momentum transport rate. The peak correlation coefficient obtained (R = 
0.25) was approximately equal to the correlation coefficient obtained Yor 
the z-0 plane and was about one-half of that obtained for the r-z plane. 

Su mmary Comments 

For this flow condition with swirl angles of 25 to 40 degrees in the 
annular stream, the principal momentum turbulent transport was in the r-z 
plane, uv, and was caused by the axial velocity gradients. The peak 
transport rates were approximately the same (0.06 m Z /s Z ) as for the 
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nonswirling flow condition but occurred nearer the inlet (z - 25 ram) than 
for the nonawu ling flow condition (z - 150 and 205 ram). The correlation 
coefficients for the momentum turbulent transport in the r-z plane were 
greater than those in the other two planes. This result is compatible with 
assumptions m the turbulent transport theory regarding the "return to 
leotropy" for eddy dissipation of energy. 


Mass Transport 

The mass turbulent transport rate profiles were obtained for the three 
measurement directions. These transport rates are generally associated 
with the gradients of the mean concentration profiles (Fig. 20). The 
discussion of the mass turbulent transport will be related to these 
concentration profiles, to the shear regions shown in Fig. 1 and to the 
flow visualization results. 

Radial Direction 

. v . The 0 l & ^\ mass talent transport rate profiles, 7f, which are shown 
ri *of * 2 ! j^ dlcate that radial transport occurred at axial locations, z = 
id, 25 and 50 mm and essentially no mass turbulent transport occurred 
between z = 102 mm -through z * 203 mm. The peak rates for z = 13 mm and 25 
mm occurred at the radial location, r/R = 0.25 which corresponds to the 
location of the interface between the inner jet and annular jet. At z = 25 
mm, the peak transport rate occurred for 0.18 < r/R <0.3 where the axial 
velocities were negative. At 51 mm, the turbulent transport rate data 
appear to be scattered about zero especially near the centerline region. 
According to the flow visualization results and the axial velocity 
profiles, the axial location of z = 51 mm is near the upstream end of the 
centerline recirculation region where the flow may not be axisymmetric and 
thus may account for the data scatter. The peak fluctuating radial. 

f ° r ewiriing flow were approximately 60 percent greater than the 
peak obtained for nonswirling flow and therefore higher mass turbulent 
transport rates should be_expected for swirling flow. However, the peak 

r /2 n r o? 8 ?° rt rate 0f Vf " °* 018 m/s which was obtained at z - 25 mm and 
r/K o 7 18 a PP r °ximtely equal to the peak value obtained for- 

nonswirling flow (Ref. 9). The peak for nonswirling flow also occurred in 
the interface region (0.20 < r/R < 0.25) but was obtained much farther 
downstream (102 mm < z < 151 mm). The peak values were relatively low 
although they occurred at locations where the radial concentration 
gradients were the highest in the flow field (see Fig. 22). For the 
swirling flow condition, the mean radial velocities were an appreciable 
fraction (20 to 70 percent) of the local axial velocity. Thus there was 
more convection of the inner jet fluid radially outward than occurred for 
the nonswirlmg flow condition. 
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The correlation coefficient, R , profiles which are presented in Fig. 
29 generally followed the trends of v the mass transport rates. The peak 
correlation coefficient occurred at the same location as the peak mass 
transport rate (z « 25 mm, r/R - 0.26) but its value of R . » 0.31 is 
approximately half the correlation coefficients obtained for the 
nonswirling flow. 

Axial Direction 


The axial mass turbulent transport rate, uf, profiles are presented in 
Fig. 30... The greatest axial mass transport (0.25 m/s < uf < 0.55 m/s) 
occurred near the upstream end of the centerline recirculation zone (r/R < 
0.3 and 25 mm < z <51 mm). In this region the mean axial velocity was 0 
negative. No appreciable axial mass transport occurred for z > 102 mm. 
There is a significant difference in the axial mass turbulent transport 
rates between the nonswirling -M swirling flows. For nonswirling flow, 
negative mass turbulent transport rates were obtained in those regions 
where the inner stream concentration was decreasing with axial location 
which indicates a counter-gradient mass transfer. For swirling flow, the 
axial mass turbulent transport rates were positive (gradient mass. 
transfer) . The peak axial mass turbulent transfer rate for the swirling 
flow (0.55 m/s) was higher than the peak obtaned for non-swirling flow 
(-0.35 m/s). The mass transport in the axial direction for swirling flow 
was essentially completed within z ■■ 102 mm whereas the mass transport for 
nonswirling flow was not completed until z - 305 mm. The axial mass turbu- 
lent transport rates were higher than the radial mass transport rates even 
though the peak radial concentration gradients were approximately seven 
times greater' than the peak axial concentration gradient. 


The correlation coefficient R f ,. profiles presented in Fig. 31 show 
the development of the axial mass transfer process more clearly where the 
mass transport rates are small. At z * 13 mm and z = 25 mm, the axial mass 
transport correlation coefficients were generally positive for r/R < 0.40 
which indicates mass transport in the downstream direction. In the shear 
region between the annular jet and the recirculation zone near the outside 
wall, 0.45 < r/R Q < 0.8, the correlation coefficient was negative for z * 

13 mm and z » 25 mm indicating mass transfer upstream. . Near the upstream 
end of the centerline recirculation zone (z < 51 mm), the correlation 
coefficient was positive near the centerline (r/R < 0.4). At the 
downstream end of the recirculation zone, (102 mm°< z < 152 mm), the 
correlation coefficient, was negative. The correlation coefficient is 
generally positive near the peripheral wall (r/R >0.8). Although the 
peak axial mass turbulent transport rate for swirling flow was greater than 
that for nonswirling flow, the correlation coefficients for swirling flow 
were less than that for nonswirling flow. 
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Azimuthal Direction 

The azimuthal mass turbulent transport rate, and the corresponding 
correlation coefficient, R -, profiles are presented in Figs. 32 and 33, 
respectively. The azimuthal mass turbulent transport rates shown in Fig. 

32 were small compared to the rates in the axial and radial direction. At 

z " ! 3 raa » ^ ere was a 8ma11 negative rate in the interface region between 
the inner and annular jet streams. This is caused by the preferential 
turbuent eddy orientation in -the interface between the inner and annular 
streams. The negative value indicates that the inner jet fluid (from the 
nonswirling inner jet) is rotating at a lower rate than the mean azimuthal 
velocity. For axial locations, z = 25 mm and z = 51 mm, there were also 
^variations of the mass transport rate near the centerline region, r/R 
< 0,4, which can be attributed to coupling between the mass and swirling 
momentum transport in those regions. For z > 102 mm, the azimuthal mass 
transport was essentially zero. 

The correlation coefficient profiles presented in Fig. 33, essentially 
showed the same trends as the mass transport profiles. There was scatter 
in the results in and near the ends of the center recirculation region due 
to asymmetries and wandering of the center of rotation. 

Summary Comments 

The mass turbulent transport process for the "simple" swirling coaxial 
ow configuration utilized in this experiment- is very complicated. From 
d *5 a fV? 6 fl0W visualization high speed motion pictures, it appeared 
that the turbulent transport and mixing occurred in several steps of axial 
and radial mass turbulent transport coupled with a large radial or mean 
convection flux. The turbulent transport rates indicated mass flux across 
streamlines and the decrease in concentration fluctutaions indicated mixing 
of the fluid toward the molecular level. These processes appear to begin 
coacentr f tlon fluid from the inner jet diffusing turbulently into 
fluid from the recirculation cell near the centerline for 25 mm < z <50 
mm. This diluted inner jet fuid (f = 0.5) is convected by the negative 
axial velocities at z - 25 mm and 0.2 < r/R <0.3 to the large eddy shear 
region at 0.3 < r/R Q <0.4 and 25 mm < z < 75 mm. There appeared to be 
more turbulent transport across the diagonal streamlines (r/R > 0.3 and 
<150 mm) by axial mass turbulent transport than by radial ma§s turbulent 
ranspor . his result may be due to the preferential turbulent structure 
developed by the shear layer. Additional analysis of. the data and 
comparison with numerical code predictions or results from dimensional 
analysis (using the equations for turbulent transport) will be required to 
fully understand this turbulent transport process. 
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DISCUSSION OF SKEWNESS AND KURTOSIS RESULTS FOR 
VELOCITY AND CONCENTRATION PROBABILITY DENSITY ..FUNCTIONS 


and transport rata ^trlWUo^"^ 1001 -’' CMlcaI '«“ i °“ distributions 
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trends indicated by the p.d.f.s are consistent with the axial velocity __ 
profiles shown in Fig. 13. 

Radial Velocity 

Although the radial velocity p.d.f.s presented in Fig. 35 showed less 
variation with radial location than the axial p.d.f.s, they did have 
several distinct features. At radial locations r/R = 0.0 and 0.13, the 
p.d.f.s were nearly symmetrical, i.e., S = 0.2 and°were sharply peaked 
near V * 0. The p.d.f. at r/R * 0.24 also peaked, about zero but had a 
negative skewness. The tails 8n each side of the peak cause the kurtosis 
to be relatively high, i.e., K 5.0. At the radial location r/R =0.36, 
the p.d.f. is skewed toward higher velocities indicating that the fluid in 
this region was being accelerated. At r/R = 0.52, the p.d.f. is 
essentially symmetric which is consistent with the radial velocity profile 
shown in Fig. 14, i.e., the mean radial velocity at this radial location is 
at the maximum. At r/R = 0.68, the radial velocity is close to zero but 
the p.d.f. is skewed to°higher velocities which indicates that some 
filaments of fluid with relatively high radial velocity occasionally passed 
through the probe volume. 

Azimuthal Velocity 

The azimuthal velocity p.d.f.s which are presented in Fig. 36 show 
that the fluid had essentially zero velocity for radial locations up to 
r/R = 0.24 (p.d.f. 8 peaked near zero). At r/R Q = 0.36, the p.d.f. has 
shifted to higher velocities. However, low velocity fluid occasionally 
passed through the probe volume since the p.d.f. is skewed toward the lower 
velocities. The p.d.f. at r/R = 0.52 where the velocity reaches a maxi- 
mum is sharply peaked and nearly symmetric (S = 0.06). At r/R = 0.68, 
the p.d.f. is skewed toward lower velocities which is consistent with the 
fact that the fluid is decelerated in this region* 

Concentration 

The concentration p.d.f.s are presented in Fig. 37 for the axial 
location z = 25 mm. Each of the profiles show unique and rapidly changing 
features which were characteristic of specific flow regions. The concen- 
tration p.d.f.s at r/R * 0.0 and 0.52 were obtained from measurements at 
the inner and outer ed§es, respectively, of the mean concentration profiles 
(see Fig. 20). At r/R =0.0, the most probable fluid concentration was 
near 1.0 with the tail°of the p.d.f. skewed toward the lower concentra- 
tions. The values of the skewness, S-, and kurtosis, K^, were -2.91 and 
12.4, respectively. At r/R = 0.52, the most probable fluid concentration 
was near 0.0 with the p.d.f? skewed toward values, up to 0.2. The skewness 
and kurtosis of the p.d.f. at this location were 3.92 and 27.0 
respectively. Note the precipitous slope of the concentration p.d.f. at 
zero concentration which is the shape ex pected for an ideal seed and 
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measurement system. The p.d.f. 's at the intermediate radial locations 
indicate the rapid mixing which is occurring in this region. At r/R ■ 
0.13, the concentration fluctuations were increasing (see Fig. 21) and the- 
mean concentration had decreased slightly (see Fig. 20). The p.d.f. at 
this location is consistent with this trend in that it is skewed to values 
of lower concentration indicating low concentration fluid was being mixed 
into the high concentration inner stream. At r/R ■ 0.24, the concentra- 
tion fluctuations were at the maximum for this axial location and the 
p.d.f. reflects this result. The p.d.f. is double peaked which indicates 
that large eddies of high and low concentration existed at this location. 
The mean concentration value was approximately 0.5 which occurs between the 
two peaks. At r/R Q * 0.36, the concentration fluctuations have decreased 
and the mean concentration has decreased further. The p.d.f. at this 
location has become more symmetric but was still skewed to higher 
concentrations which indicates that the mean activity had decreased but 
high concentration fluid was occasionally being carried into the low 
concentration region. The p.d.f. at r/R ■ 0.68 is essentially symmetric 
and sharply peaked at a value of approximately 0.05 indicating very little 
concentration .fluctuation. The p.d.f .s at higher values of r/R are not 
shown but were also sharply peaked with the peak near a value o? 0.05, the 
mass flow averaged concentration level within the annular recirculation 
cell. 


Skewness and Kurtosis Distributions 

The skewness and kurtosis distributions for the axial, radial and azi- 
muthal velocity p.d.f.s and the mean concentration p.d.f. are presented in 
Figs. 38 through 45. The skewness, S, is a dimensionless measure of the 
asymmetry of a p.d.f. If the p.d.f. is symmetric about the origin then the 
skewness would be zero. A p.d.f. with positive skewness indicates that 
there are more values of the function in the tail toward the positive side 
while the negative skewness indicates that there are more values of the 
function in the tail toward the negative side. The kurtosis, K, is a 
dimensionless measure of the flatness of the p.d.f. The value of kurtosis 
is large if the values of the function in the tails of the p.d.f. are rela- 
tively large. The kurtosis of a Gaussian p.d.f. profile would be equal to 
3.0 while the kurtosis for a square p.d.f. profile would be 1.8. 

As will be noted in the discussion which follows, the skewness of the 
velocity p.d.f.s generally fell within a band of values equal to +. 2. The 
kurtosis factors were generally greater than 3.0. On the other hand, skew- 
ness factors for the concentration p.d.f.s were as high as + 10 with 
kurtosis factors approaching 100 in some flow regions. It has b een sh own 
(Ref. 15) that utilization of Schwarz's inequality, i.e.,(x?) 2 < x 4 x 2 with 
the definitions for skewness and kurtosis (see Appendix III) leads to the 
relation 
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K > S" 


(3) 


fers it 8 f:„T;: s j t 5 i V l r 8e A8 in B r?f n i! ude <of ai ; h « ai * a >- *• 

£ono»s. the kurtosis KlSS; tMe^e^ali"? “ the »*ich 

Axial Velocity 

ted in Fig. 38. For°axial’ locations l P *d* f * 8 *** pre8en " 
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to scatter about zero Whenever n,. u ame UI ) lfor ®> the skewness tended 
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tion (Fig. 39). * * e vadue ^ or a Gaussian distribu- 
te s 2: *■«««. .th. t 

velef : t appea ^f t0 be C0rreUt6d t0 te local curvature of the Lui ,ai “ 
elocity profile. It appeared that. S u > 0 whenever 3 5 U/ 9r 2 > 0 and S < 0 

whenever 3 U/3r < 0. These results have been explained (Ref isu” 

^ocal 8 ^ ^ St a radial location where the mean velocity profile has^ 
local maximum, slower movine flniH ui'ii ^ • ,, ^ prorne has a 
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Radial Velocity 

The skewness and kurtosis distributions for the radial velocitv 
Ped.f.s are presented in Fies 4f) anti At >ru , auiai velocity 

skewness occurred at the Lial w a H n ! "i The largest variation in 
section, s - 5 » F« i^eLna ^l *?* entra "“ of tha «« 

became flatter and the Knde! to ... '.“‘“i the SkaM " eas profilaa 
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upstream end of the centerline recirculation region? r/K, < 0 ” 

Azimuthal Velocity 

Although the .skewness distribution, for the azimuthal velocity p.d.f ., 
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shown in Fig. 42 indicated that there was considerable scatter of the data, 
the average values of the skewness factor, S , are generally compatible 
with the local fluid acceleration or deceleration. The kurtosis profile 
shown in Fig. 43 also exhibited scatter in the data but except for some 
isolated radial locations near the entrance of the test section, z ■ 5 mm, 
and well downstream of the entrance z > 305 mm, the kurtosis values 
generally averaged between 3 .0 and 4.0. 

Mean Concentration 

The skewness and flatness factor profiles obtained for the mean 
concentration p.d.f.s are presented in Fig. 44 and 45 and behave as 
expected in most of the flow regions. In the high concentration region (z 
=13 mm and 0.1 < r/R < 0.25), the skewness, S f , was negative indicating 
that eddies of low concentration fluid are occasionally passing through the 
probe volume. For those regions where the concentration gradients were 
large (z » 13 and 25 mm; 0.2 < r/R <0.4), the skewness rapidly changed 
from large negative values to large positive values. Even in the mixed 
regions, (r/R > 4; f 88 0.1), the skewness was greater than zero because 
some filaments of fluid with relatively high concentration occasionally 
passed through the probe volume. Farther downstream in the test section, z 
> 152 mm, the skewness leveled out to essentially zero. 

The kurtosis profiles presented in Fig. 45 exhibit similar trends. 

Very large values of kurtosis, K f > 100, occurred in those regions. where 
the concentration gradients are steepest. As the fluid mixing was 
completed, the kurtosis profile flattened out to an average value of about 
4JL 
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Typical Turbulent Transport Rate Probability Density Functions 

transport P rate b datrarrpJeLnted t forthe r a^ialT e "t“'“ “ d ” a88 tur '“ ulent 

IS the same location for which the velocity and J° Catl ° n * * " 25 mm » which 
presented* This location W was C c2l« # i; C 52iSl"S n Jh nt « ti0n p : d,f ' 8 were 
the mass and momentum turbulent transfer rates arf the highest? 1011 WherC 

Momentu m Transnm-t- 

in the^r-z, f ° r ^ ® omentum transport rates 

and 48 respectively. The d d f o C pla ? e ® are . presented in Figs. 46, 47, 

distributions about the zero momentum turb^wV* 1 ** they * 11 had peaked 
typical of the shape of momentum rat! « a transport rate which is 
(Ref. 16).. The p.d.f.s for 1 r-6 olL turbulent boundary layers 

those for the r-z and z-0 planes A *! Were m0re sharpl y peaked than 
essentially the same variation in ° f profiles showed 

changes in shape occurred at those radia^loelt • ^ 1 ? cation * The greatest 
turbulent transport rates were the h !! f Rations where the momentum 
0.68. 8 W6rt the kl 8 he st j l .e. , r/R = 0.24, 0.36 and 


Mass Transport 


transport P rates which^J^pCefStCd 0 ^ Fig^/*^* 1 mass turbulent 
profiles are sharply peaked near zer! ^ lnd ^ cate that all of the 

< 0.52 are broader indicating incr^d ^ill*’ P ‘ d - f - 8 3t 0al3 < 
this region. This trend is Consistent wiJhfh, maS8 transport activity in 
presented in Fig. 28 which show that thl ?V h ? m38S trans P° rt rate result 
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Probability density functions of the axial maBB „ . i 
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Typical Momentum Turbulent Transport Results 
Radial^Axial Plane (r-z) 

The second central moment profiles (or rms fluctuation from the mean) 
of the momentum turbulent transport in the r-z plane, a , are presented in 
Fig. 52. This quantity previously was used to analyze aXd evaluate the _ 
turbulent transport process in boundary layers, e.g., Ref. 17. A 
comparison of the values of a with the corresponding values of the 
momentum transport rate (57) ffom Fig. 22 show that a was generally at 
least a factor of two greater than uv. Ratios of 7uv equal to 
approximately 3 were previously reported for boundary layers (Ref. 16). 

The trends shown in Fig. 52 were generally the same as shown in Fig. 22. 
highe^X W6re Senerally highest were the momentum transport rates, uv* were 


Skewness and kurtosis profiles (S and K ) for the momentum turbu- 
turbulent transport p.d.f.s obtained in the r-z v plane are presented in Fig. 
53 and Fig. 54 respectively. The skewness profiles showed the same trends 
as the momentum turbulent transport rate. The skewness was positive when 
the momentum transport rate profiles was positive and became negative when 
the momentum transport rate became negative. For those axial locations 
where the momentum transport rate became zero, the skewness values tended 
to scatter about zero. 

The kurtosis profiles presented in Fig. 54 showed the trend expected 
from the skewness profiles. Whenever the skewness was large in magnitude 
of either sign, the kurtosis was large and positive. However, the low 
kurtosis values obtained at z = 5 mm, r/R = 0.45 indicate deviations from 
the processes assumed in the simple transport models such as the. joint- 
Gaussian p.d.f, model. J 

Axial Azimuthal Plane (z-r) 

• ... The ® ec ° nd central moment profiles of the momentum turbulent transport 

m the z-0 plane, cr , are presented in Fig. 55. The .o profile generally 
peaked where the momentum transport rates (17) peaked bu£ were 2 to 5 times 
greater in value. As the momentum transport rate dampened to zero the 
T® 1 "® 8 of °uw also approached zero. The skewness profiles of the p.d.f.s 
the z~6 plane which are presented in_Figs. 56 generally followed the 
same trends that were observed for the uw profiles. The skewness was 
generally positive when the turbulent transport rate was negative. 

owever, many of the fluctuations which occurred in the momentum rate were 
not observed in the skewness profiles. The kurtosis profiles, K 
presented m Fig. 57 exhibit a large amount of -scatter.. However? w the 
kurtosis tended to be large whenever the skewness was large in magnitude. 


35 


R83 -91 5540=26 


Rfldial -Az imuthal Planp (y-fll 
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Summ ary Comments 
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Typical Mass Turbulent Transport Results 
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profiles. At axial locations z * n mm 8 turbulent transport rate 

central moments occurred in the interfL* J z .“ nBn ' the Peak second 
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th-n f f adla i ma88 turbulent transfer rates were slightly higher 

than the peak values for the azimuthal mass transfer rates. The peak 
second central moment for the axial mass transport is approximately 1 5 
times the peak observed for the radial and azimuthal mass transfer! Also 
the peak values of the second central moments for all measurement 

? enerall y higher than the peak mass transport rate for any 
direction at that measurement location* ' 

p . 1r ne r^ f * the u S8S turbulent transport rate p.d.f.s are shown in 

Figs. 62, 65 and 08 for the radial, axial and azimuthal directions 

respectively. The skewness values deviate from zero where mass transport 
occurs but no trend consistent with the mass transport p^ofUes is P 
apparent • 

Fie 6?%^ r ^ 8i8 of . the ® ass transport p.d.f.s are presented in 

Fie* 69 for T l & * J d 1 lr f tlon ! in Fi 8* 66 *°r the axial direction and in 
rh?** 69 5 the azimuthal direction. These kurtosis profiles also followed 
th. trend, of the kurto.is profiles discussed previously. Whenever the 

b^n.rgr e Ur86 “ " a8 " itUde ( ° £ eUh » *•>. the y kur tostn Uo 
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SUMMARY OF RESULTS 


Qualitative and quantitative studies were conducted of the flow down- 
stream of swirling coaxial jets discharging into an expanded duct. The 
ratio of annular jet diameter and duct diameter to the inner jet diameter 
were approximately 2 and 4, respectively. The inner jet peak axial 
velocity was approximately one-half the annular jet peak axial velocity and 
the mean swirl angle in the annular stream was approximately 30 degrees. 
Results from the studies were related to the five shear regions within the 
duct: (1) the wake region downstream of the inlet, (2) the shear layer be- 
tween the jets, (3) the annular recirculation region, (4) the reattachment 
region, and (5) the .centerline recirculation region. 

A flow visualization study was conducted using fluorescence dye as a 
trace material and high-speed motion pictures to record the dye patterns in 
selected r-z and r- planes. The results of the flow visualization study, 
are summarized with the following observations: 


1. The flow was as axisymmetric as could be determined visually... 

2. The larger scales of the turbulent structure were observed to grow 
within the centerline recirculation region. This growth occurred from- 
the width of the wake region downstream of the inner jet tube to a 
large fraction of the duct diameter. 

3. Downstream of the centerline recirculation region, the flow initially 
had large lobes which evolved into a vortex swirl pattern. There was 
little radial mixing downstream of the recirculation cell. 

4. The high intensity turbulent eddies in the shear layers were not 
axisymmetric or periodic. The large scale waves and eddies appeared to 
have a range of wavelengths. 

5. Two major mixing regions were observed: (1) at the interface between 

the inner jet and the centerline recirculation zone, and (2) at the 
interface between the inner jet and annular jet streams. 

6. Mixing at the interface of the inner stream and the recirculation cell 
diluted the inner jet concentration and the resulting mixture of fluid 
from the recirculation cell and the inner jet was entrained by the 
swirling annular stream. 


A detailed map of the velocity, concentration, mass turbulent 
transport rate and momentum turbulent transport rate distributions within 
the test -section was obtained to provide data for the evaluation and 
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improvement of turbulent transport models r i * 

win.Tr ° btain H d si " ul “"<=»“= IV deumine the momLtu^tirbuUnr”"”' 

transport rate and mean velocities. Data sets of ^ urDui ent 

»aB.Mrt S r''t re < ’ bCained ““ultaneously to determine mass turbulent 6 ”*™" 

Following are the principal result, from thie study. 


The axial, radial and azimuthal velocity profiles described the 
changes m .the shear -regions within the duct. 

V™ fluctuating concentration profiles described the inner 

jet fluid distribution within the duct. 

Mixing for swirling flow was completed in one-third the length 
required for nonswirling flow. length 

The momentum turbulent transport rate measurements in thp r . ,» . 

r-9 p a„ e . da scribed tha local momentum fuSTS 0 S bu".nt ““ 

“S. a:d rr datr:« Me££iCien “ " erc ° b£ained £or measurement 

11. The principal momentum turbulent transport was in the r-* n i a „ . 

turbulfnt a transport ted t t0 th * a * ial Vel ° city gradients. Peak momentum 

sssisL, ss 0 ap ™ teiy «• - £0 r the 


7. 

8 . 

9. 

10 . 


SdS^ar^ci^ 1 ^ ^rL°nng i8 ffo r : dien xL r :: h k r tha ? counter - 
*2 srvr 8rea f er than the ^ - » t ma : t 8 es tran8 - 

ZlnllT ^ia a r 8 a ^ur ati0n gradientS — approximately^ne- 

The mass turbulent transport process for swirlin* m.viai 
complicated. Mining appears to occur il .“^.tep^f 
radia! mass transport coupled with a large radial “£n invective 

^L t rte‘nnEeJu^rrecin a utUn b a!ne “^LhdUut'T™"" fl “ id 

is then convected b, the negative Eui veUcitle * L^Ura ““ 
eddy shear region between the inner and annular .streams . 


13. 


14. 
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peak mass transport correlation coefficients obtained for nonswirling 
flow although the axial mass transport rate for swirling flow was 
greater than that for nonswirling flow. 

16. The skewness of the axial velocity, p.d.f.s was related to the 

curvature of the axial velocity profiles. S < 0 was obtained for 
3 U/^r <0 while S > 0 was obtained for 3 / U/8r%> 0. The skewness 
was also proportional to the magnitude of 


!/»* 2 ' 
9 U/3r , 
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APPENDIX I 

FLOW VISUALIZATION FOR ALL FLOW CONDITIONS 


Flow visualization studies were conducted prior to the selection of 
the flow condition for detailed data acquistion to determine the effects of 
the velocity ratio, U./U , on the flow characteristics within the test 
section. 

Motion pictures were obtained in the r-z plane with the center of 
illumination at z = 100 mm and in the r-0 plane at selected axial 
locations. Motion pictures with a frame speed of 500 per second were 
obtained for each of the following five flow conditions: 


Flow 

Mean Axial Velocity, m/s 

Flow 

Rate, gpm 

Condition 

Inner Jet, 

Annular Jet, U 

a 

Inner Jet 

Annular Jet 

1 

0.52 

1.66 

6.2 

52.8 

2 

0.27 

1.66 

3.2 

52.8 

3 

2.08 

1.66 

24.6 

52.8 

4 

0.94 

1.51 

11.1 

48.0 

5 

0.94 

2.87 

11.1 

94.8 


In the following paragraphs, the flow visualization photographs for 
Flow Conditions 2 through 5 presented in Figs. 70 through 73 will be 
discussed. The discussion of Flow Condition 1 (Fig. 10) will be repeated 
for reference. Each figure contains two photographs taken in the r-z 
plane. The photograph at the top left side shows the flow characteristics 
when dye is injected into the inner stream. The photograph at the top 
right side is a visualization of the flow field when dye is injected into 
the annular stream. The four photographs presented at the bottom of each 
figure show the flow characteristics in. the r-0 plane at selected axial 
locations when dye is injected into the inner stream. 


Flow Condition 1 

Photographs taken for Flow Condition 1 in both the r-z plane and at 
selected axial locations in the r-0 plane are presented in Fig. 10. The 
photographs taken when dye is injected into the inner stream (top left), 
shows that a high inner jet stream concentration persists for an axial dis- 
tance of approximately 70 mm.. In the photograph taken when dye is injected 
into the annular stream (top right), dye has been entrained near the 
centerline between an axial distance of approximately 60 mm to 170 mm indi- 
cating the presence of a recirculation Cell near the centerline. The exis- 
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tence of this recirculation cell becomes more apparent in the high-speed 
films and appears to extend to almost 150 mm. Both photographs show large 
scale eddy structure in the shear layer between the jets starting at an 
axial distance of approximately 30 mm and extending downstream to at least 
75 mm. The scale of these eddies at z ■ 50 mm is more than half the inner 
jet diameter of 30 mm. Although the eddy sizes are large, the eddy struc- 
ture did not appear to be periodic or symmetric. When dye is injected into 
the annular stream, the annular recirculation zone adjacent to the test 
section inlet plane and the attachment region downstream of this annular 
recirculation zone are clearly defined. The length of this annular recir- 
culation cell is approximately 50 mm which is about one-third the length 
observed with nonswirling flow for the same flow condition. 

The photographs taken in the r-0 plane at z = 25 mm shows the radial 
scale ( = 6 mm) of the eddy structure at the upstream end of the shear 
region. At z = 51 mm, the scale of the large eddy structure has increased 

818 ™n C<mtly but does not a PP ear t0 change appreciably from z = 100 mm to 
z = 200 mm. As shown by the photographs in both the r-z and r-0 planes, 
edd y structure is not axisymmetric or azimuthally periodic for 


Flow Conditions 

The inner jet velocity,. U for Flow Condition 2 was approximately 

that ~ 0W Condltion 1 * The Photographs, (Fig. 70) indicate 

that the overall flow characteristics are similar to those observed for 
Flow Condition 1. However, as a result of the reduction in inner jet 
velocity, the distance over which the high inner jet fluid concentration 
extends from the test section inlet has been decreased from approximately 
/U mm to 35 mm* J 


Flow Condition 3 

^e mean inner jet axial velocity for this flow condition was approxi- 
mately four times that for Flow Condition 1. As a result, the inner jet 
stream was flowing faster than the annular stream. The ratio of the peak 
inner jet to the peak annular velocity was estimated to be approximately 

it F l ?*-. 71 ’J the flow field for thiB “ow condition differs 
from Flow Condition 1 m that the flow becomes asymmetric and the inner jet 

f?nfnrTh a T * 6 cent ?fl ine for ^ > 70 mm. The flow appears to be in- 

enced by two forces: (1) the tendency of the inner jet to flow 

downstream without recirculating and (2) the strong tendency of the swirl- 
ing annular flow stream to form a centerline recirculation cell. The 
result appears to be a combination of the two effects in that the inner 
stream is defiected by the swirling flow in the centerline recirculation 
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Flow Condition 4 

This flow condition has approximately the same total flow as Flow 
Condition 1, however, the ratio of the inner stream mean velocity to 
annular stream mean velocity has been increased from 0.32 to 0„64. Based 
on measurements in Ref. 8, the peak axial velocities in the inner and 
annular stream are expected to be approximately equal. The photographs 
(Fig. 72) indicate that the flow field appears to be axisymmetric and 
stationary. Inner jet velocities 20 to 30 percent greater than that for 
Flow Condition 4 are required to obtain the precessing center jet 
characteristics of Flow Condition 3. Compared to Flow Condition 1, the 
length of the annular recirculation cell is 30 to 50 percent greater and 
the distance from the inlet plane to the inner recirculation cell is 
increased by 50 percent. The large scale eddy characteristics also 
appears approximately the same for Flow Conditions 1 and 4. 


Flow Condition 5 

The velocity ratios for Flow Conditions 5 and 1 are equal: however, 
the absolute velocity for Flow Condition 5 is increased by a factor of 1.8 
compared to Flow Condition 1. The photographs presented in Fig. 73 indi- 
cate that the flow characteristics for Flow Condition 5 were essentially 
the same as those observed with Flow Condition 1. 

The flow visualization study showed (1) an additional flow region: 
namely, a centerline recirculation region,, which occurs when swirl is 
introduced into the annular stream of confined coaxial jets and (2) the 
inner jet will precess when the ratio of inner jet velocity to annular 
stream velocity exceeds a critical value. These characteristics of 
confined swirling jets have been previously recognized in the literature. 
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APPENDIX II 

SWIRLER EXIT MEASUREMENTS 


Laser velocimeter measurement techniques were used to determine the 
mean and fluctuating velocity distributions approximately 5 mm downstream 
of the swirler exit face when the trailing edge of the swirler was reposi- 
tioned at the test section inlet. Data was. obtained at various azimuthal 
and radial locations to define the velocity distribution at the swirler 
exit. This survey was made because the flow field downstream of the swirl 
vanes will contain wakes caused by the blade skin friction. Also, the flow 
field will be influenced by secondary flows which result from the 
interaction of boundary layers on the vane end walls with the pressure 
gradient in the vane passage. 

Prior to the detailed data acquisition, axial and azimuthal velocity 
data was obtained at the swirler location r/R = 0.40 by rotating the 
swirler and the center duct in 15 deg increments while maintaining the 
laser velocimeter probe volume at the same location. The azimuthal 
velocity profiles were to be used to determine the swirler asymmetries and 
to select an azimuthal location for a more detailed velocity survey. Since 
the swirler has eight blades, three measurements were obtained downstream 
of each blade passage. Although the velocity distributions obtained showed 
some blade-to-blade variation, the variations were only a fraction of the 
variation across each blade. The. blade passage chosen for the detailed 
measurements was one with a local velocity distribution representative of 
the entire swirler and was the center passage of three passages which had 
nearly identical velocity distributions. This blade passage was also the 
passage in which the detailed measurements presented in the main body of 
this report were taken. 

The flow field 5 mm downstream of the selected vane passage was 
characterized by making a velocity survey of the axial, radial and 
azimuthal velocities along radii starting at r/R = 0.225 (the location of 
the thin metal tube separating the inner and annular jet streams) and 
ending at a radial location r/R Q =0.6 which is approximately 7 mm outside 
the annular passage. Data was obtained at ten azimuthal locations, five 
(5) deg apart by rotating the center tube and swirler as shown in Fig.. 74. 
The striped blade shown in Fig. 74 is a reference blade chosen for swirler 
orientation purposes. The detailed measurements reported in the main body 
of the report were taken with the leading edge of the reference blade set 
at <|) = 75 deg. The LV measurements were made for. the axial/azimuthal 
velocities by traversing horizontally first, along the radius indicated by 
swirler orientation <J» * 50 deg (leading edge of reference blade set at <f> = 
50 deg). Horizontal traverses were then made along the other radii indica- 
ted in Fig.. 74 by rotating the reference blade through 5 deg increments. 

In this way, the flow field downstream of a single swirler vane, and vane 
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n*ffw e ^? U i d b f cba ? acterized * Lv measurements were made to determine the 
axial /radial velocities by first rotating the center tube and swirl™ 90 
degr and then traversing vertically alone each radii Th* -* r 

th! trailing LgTof po^SedT** 

iith the trailing ad« (R “ n ' 59 throu8h 78>i snd then repeated 

Tnioi. i trai “ n 8 ed 8e of the swxrler vane positioned 51 mm upstream of the 
inlet -plane of the test section (Runs 79 through 98). P the 

tations f-lo^S^O^nd^ r l0City mea8u ? ements made a * airier orien- 
caciona <(. - 50, 55, 60 and 95 deg were made in the passage between blades 

effe«s"o? M.d th ? “ Ul Profile no 

art at bl d ? k ??‘ f 180 ’ 8ince 8wirler orientations 4 = 50 and 95 deg 
are at the same relative location for two vanes, it would be expected that 

identiial lty Th 18tr ^ b ? t10 ? 8 -° r the8e two locations would be nearly 
these expectations^^ ’ y Pr ° fiUS > >reBented in Fi 8- ” substantiate 

i c« 8 8 j°^ j 11 the velocity distributions at swirler orientation 

* 50 and 95 de 8 are nearly identical and are easily identifiable when 

compared to the intermediate orientations, as being periodic profiles The 
l?*} 1 ” 0b 5 alnedat * " 5 5 and 60 deg are similar and shoilittleeffect 
en bl de tr ke8 *,.?° wev ! r ’. the effect of b lade wake on the axial velocity is 
dee Th^ he Pr0fl J eS obbained for swirler orientations <j) = 65 through 90 
deg. The vane wake produced a dip in the axial velocity profile which 

™ d ?* 3 dif 5 6rent radial ^cation for each Airier Irlentltlo* Note 
c^e s^n 8t IP occurred for a swirler orientation of 80 deg which 

bU^Uee Fig. m 7 e 4r re0ent8 e88entiall y a l™8 the central part of the 


v . f an radial velocit y Profiles are shown in Fig. 76. The radial 

velocity downstream of the inner jet increased from a value of zer^aJ the 

6 b ?!‘ llne ’ as ex Pected for axisymmetric flow. Small differences in the 
profiles occurred for 0.22 < r /R < 0 .25 and can be attHbu^ to 

di^^w^ ° f c th u inne 5’ fl °" with the flow from the swirler. The 
istribution of the radial flow inward and outward for 0 25 < r/R < n s • 

attributed to (1) the formation of corner vortices due to the Inv?sci2' * 

vanes r and X (2)°the°radiall 88 ** fl " thr0U8h the ^ the 

blSI wakesf For S tT:”^^” l\ X Zll *1 ^ 

pr ° t aUs essentially smooth^and^free^of ^th^ef feet s of 
S secoUarfv^tJce's"^ 81 Vel °' Uy Pr °' iU8 « F ~> 8 ^ offered' 

The azimuthal velocity profiles presented in Fig. 77 show lareer 
relative velocity changes, due to the'blade wakes a^d secondary ?Iows in 
the vane passages, than the axial and radial velocity orofil*/ Th ■ 
tions in the azimuthal velocity profiles due to tK hl5 * l J Vari8 ' 

the same combinations of suirl^££ti“ the":"^ “ 
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velocity profile, i.e., $ * 60 deg and r/R„ - 0.3 to <j> » 85 deg and r/R « 
0.45. Note the azimuthal velocity was slightly greater for 0.5 < r/R $ 

0.6 than for r/R Q « 0.48. In this region, 0.5 < r/R < 0.6, the axia? 
velocities were near zero and the flow was coming from the recirculation 
£ el l * 18 shows that the angular momentum imparted to the flow 

V' R o ' was ^ ess ^at in the outer recirculation cell 

and hence less than the average angular momentum imparted over the entire 
blade. This may be due to the leakage around the outer tip of the swirler 
blades. The blades are soldered to a hub which is fastened to the inner 
annular wall but they are unattached to the outer annular wall and have 
approximately 0.1 mm clearance. Some variations in the azimuthal velocity 
profile, e.g., at d> » 75 deg and r/R ■ 0.4, may be caused by the secondary 
flows developed in the vane passages u due to the axial velocity profile at 
the vane entrance. Another interesting feature of the azimuthal velocity 
profile is the negative velocity at <t> ■ - 80 deg and r/R « 0.47. The 
azimuthal velocity profiles also show variations at raSial locations 
where the axial velocity profiles were smooth. With an almost uniform 
axial velocity profile in the center of the vane passage, the radial and 
azimuthal velocity profiles should be coupled by any secondary vortex flow. 

The axial, radial and azimuthal velocity profiles obtained with the 
trailing edge of the swirler vane was positioned 51 mm upstream of the test 
section inlet plane are presented in Figs. 78, 79 and 80, respectively. 

The axial velocity results presented in Fig. 78 show that the profiles were 
more nearly uniform but the effects of the blade wakes were still present 
although the dips in the velocity profiles associated with the wakes -were 
not as pronounced. The radial and azimuthal velocity profiles shown in 
Figs. 79 and 80 also show the same trend. It appears that even though the 
swirler is set back 51 mm from the inlet plane of the test section, the 
velocity profiles a short distance downstream from the inlet to the test 
section were not axisymmetric because the effects of blade wakes and/or 
secondary flows still persisted. These asymmetries in the flow were 
apparent in the velocity profiles obtained at an axial location of 5- mm 
which were presented in Figs. 13, 14 and 15. 


The fluctuating axial, radial and azimuthal velocity profiles obtained 
when f.ne trailing edge, of the swirler was positioned at the inlet plane are 
shown in Figs. 81, 82 and 83, respectively. Note that the turbulence 
intensity in the core region was approximately 3 to 4 percent of the axial 
velocity which is the magnitude expected in the central region of developed 
tube flow. The fluctuations were highest at the interface between the 
inner and annular jet (r/R ■ 0.25) and between the annular jet and the 
recirculation zone near the outside wall (r/R = 0.5). There were also 
relatively high fluctuations in the blade wake regions. 

. lu 6 f lu ?5V at * n8 axial, radial and azimuthal velocity profiles obtained 
when the trailing edge of the swirler was positioned 51 mm upstream of the 
inlet plane of the test section are presented in Figs. 84, 85 and 86 
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respectively. The fluctuations at the interfaces were still present but 

reduced e.p.ci.lljr et the interfece betyeeS 
fi not - Jet , the an ? ular J et * The greatest reductions occurred in the 
tuation^due It Substantial reductions in the flue- 

c“mpi“ely dr“«p'pe.“ J>CCUrred but the ““ »°t 

The momentum turbulent transport rate distributions in the r-z plane 
ini * hl f h were obtained for the swirler exit plane at the test section 
inlet plane are presented in Fig. 87. This transport is related to the 
“"l" '«•' W 8r ,‘ For those profiles without blade wakes Jeer 

5 50? ,1 ”" lrl ! r . (r/ ?o.- °' 25 > « the OD of the swirler (r/R . 

0.507, the transport is similar to that shown in Fig. 22. For the swirler 
orientations where the blade wakes are present, the momentum turbulent 

Note^hlt W h r ? ° f ^ 8ame magnitude as th08e in the ID or OD region. 

Note that the turbulent transport rates in the core regions were 

negligible. The momentum turbulent transport rate distribution, 

o tamed when the swirler was moved 51 mm upstream of the entrance to the 

test section are shown m Fig. 88.. The profiles have become more uniform 

and the momentum transport has been substantially reduced in the blade wake 

£6§X0IIS«_ 

m Fif e 8r;^r t ^ bule r- tran ? por ^ rates in the z_e pune ’ are 3h ° wR 

M/ll Thl It ! ran ! port 13 related to the strain rate: (l/r)3U/36 . 

o W/ 3 z . The momentum transport was highest at the interface between the 

iake r revion nnUl ?^-' ,et8 8i8nificant transport also occurred in the blade 
factor 2 to ™ 18 momaR tum transport rate was substantially reduced (by a 
factor or 2 to 3) when the swirler was moved 51 mm upstream of the test * 
section entrance as shown by the distributions presented in Fig 90. The 

P I° ll6b * n . F * 8, 9 ? are veR F similar to the profiles presented in Fig. 24 
at the axial location z « 5 mm. 8 H 

The profiles presented in Figs. 75 through 90 show that significant 
hree-dimensiama variations can occur in the flow field immediately 
downstream of a typical combustor like swirler. They also provide insight 
to initial conditions required for the three-dimensional swirling flows 
downstream of a swirler inlet. The data was obtained in a slmultaneoustwo 
velocity component data acquisition mode in the same manner as the turbu- 

The^f or! d ?- 8 Wh r Ch WaS pre8ented in the ma i« body of this report. 

of f ^ % h6 i flr8 ° Ur mon,ent8 of the three velocity components as well 
It- £ a t* -° Cal momentu ® transport in the r-z and z-e planes were obtained. 
jJ 1 - data 18 Presented m tabular fora in Tables IV-59 through Tables 
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ORIGINAL PAQIHS 
OF POOR QUALITY 


APPENDIX III 

DEFINITIONS OF SKEWNESS AND KURTOSIS FOR 
VELOCITY, CONCENTRATION, AND TRANSPORT 
PROBABILITY .DENSITY FUNCTIONS 

Terms in this appendix for the velocity components and concentrations are. c'._.inc'd 
using the notation of Ref. 18 and conventional statistical methods. 


B(u) 


o u or u 




Ku 


Local instantaneous axial velocity component 

Probability density function (p.d.f.) of u with properties B(u) > 0 
and B(u) du B 1,0 

Mean value of axial velocity component defined: U = u B(u) du 

Local instantaneous axial velocity fluctuation from the mean, defined: 
u.= u - U 

Second central moment of velocity u defined: a ^ = u'^ = / -c ” u^ B(u) du 

U —CO 

Will also be denoted as rms fluctuation. 

nth central moment of velocity_u defined: u n = u n B(u) da 

Skewness of velocity component, u,. p.d.f. defined:.- S u = u^/a u ^ 

Kurtosis (or flatness factor) of velocity component, u, p.d.f. defined: 


uV o 


u 


In like manner, the mean, rms fluctuation, skewness, and kurtosis for the radial 
velocity, azimuthal velocity and concentration are defined.. 


The second moments, skewness and kurtosis for the momentum and mass transport 

rates are defined in a similar manner. 

uv- ■ Local instantaneous momentum turbulent transport rate: Cu-UHv-V) 

B(uv) Probability density function (p.d.f.) of uv with properties B(uv) > 0 

and /l 00 B(uv) d(uv) ~ 1.0 

Ou 

uv Mean value of turbulent momentum transport rate defined: uv = /1’ 0 > (u-U) 

(v-V) b(uv) d(uv) 

(uv)' Local instantaneous fluctuation of momentum transport rate from mean, 

defined: (uv)' = uv - uv 
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°uv Second central moment of momentum -transport rate; 

^uv " /"„-(yvV Z B(uvj d(uv) - 

(uv) n nth central moment of momentum transport rate:. 

(uv) n * (uv) ,n B(uv) d(uv) 


uv 


K, 


uv 


Skewness of momentum transport rate: 


Kurtosis of momentum transport rate: 


s uv * <iv) 3 /^ v 

K uv - (Gv) 4 /o 4 • 
uv 


In.a 

momentum 

defined. 


like manner, the mean, second central moment, skewness and kurtosis 
transport m the r-z plane and the mass transport in three directions 


for the 
are 
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TABLE I 

Components Used In LV and LV/LIF Measurement System 

I. Laser Light Source 

Argon Ion Laser (Spectra Physics Model 164) 

TEM„ mode 

All lines, 1.0 watt power 

0.4880 urn wavelength, 0.5 watts power. 

II. LV Optics 

DISA Type 55x00 Two Color LDV System 
Polarization rotator 
Beamsplitter - Module. I 
Bragg Cell, 1 mHz effectiv£ frequency offset 
Beamsplitter * Module II 

Backscatter Section - 0.4880 pm wavelength filter and photomultiplier tube 

Backscatter Section - 0.5145 pm wavelength. fi.lt.er ancLphotomul tipi igr tube 
Pinhole Section ^ 

Beam Translator 
Beam Expander 

Achromatic lens, 310 mm focal length 

III. Electronics 

LV Signal Processor (SCIMETRICS Model 800A) 

2 units 

0.4 to 2.0 mHz range 
3% data window 

4/8 and 5/8 comparison for "good signals" 

Oscilloscope (Tektronics Model 465B) 

2 units 

LV Data Handling Interface (UTRC design) 

Clock 

Coincidence check 
Minicomputer (DEC POP 11/10) 

Floppy disk 

DECwritir III (1200 baud rate) 

IV. LIF Electronics 

Low Pass Filter (Kronhite Model 3202) 

2 KHZ 

Voltage Amplifier (Preston 8300 XWB Amplifier - Model A) 

1-1000X Amplification 
A/D Converter (DEC LPS11) 

Computer Controlled 

Digital Voltmeter (Hewlett Packard Model 3465A) 

High Voltage Power Supply (Fluke Model 415B) 

0-2500 volts 
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Test Date: 10/7/82 Run No.: 28 Flow Condition: 1 Geometry: 

Axial Location: 25 mm (1.0 in.) x/R = 0.416 
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AXIAL AND RADIAL VELOCITY DATA AND CORRELATIONS 
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FIG. 


R 83-9 15540-26 


ORIGINAL PACSti ft] 

OF POOR QUALITY 


SHEAR REGIONS WITH CONFINED, EXPANDING COAXIAL JETS 


a) SWIRLING FLOW 

LARGE EDDY SHEAR - 
WAKE REGION BETWEEN JETS 



b) NONSWIRLING FLOW 

WAKE 



RECIRCULATION 
REGION ... 


FULLY 
DEVELOPED 
DUCT FLOW 


FULLY 

DEVELOPED 

DUCT 

FLOW 


83 — 7 — 19 — 1 " 



SCHEMATIC OF TEST FACILITY 














OPTICAL ARRANGEMENT FOR FLOW VISUALIZATION PHOTOGRAPHS AND MOTION PICTURES 

IN r-z PLANE 




PICTURE CAMERA 


R83-91 5540- 26 


ORKSffJAt Pfi'ZC SST 

OF. POOR QUALITY 
















RECEIVER BEAM PATH 












VISUALIZATION OF FLOW CONDITION 1 WITH SWIRL IN ANNULAR STREAM 



bYE ADDED TO INNER STREAM 





883 - 915540^-26 


ORIGINAL PA0E Kfl 
OF POOR QUALITY 


VISUALIZATION OF FLOW CONDITION 1 WITH SWIRL IN ANNULAR STREAM 


U, s= 0.52 m/s U a = 1.66 m/s 
Q| = 6.2 gpm Q a = 52.8gpm 

DYE ADDED TO ANNULAR STREAM 

(■8 PLANE 


z = 25 mm 


z = 38 mm 






2 — 51 mm 


z = 102 mm 


2 = 152 mm 




z = 203 mm 

z ~ 305 mm 

z - 406 mm 



83 - 7 - 19-31 


a 
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R83-*5l 5540*26- 


ORIGINAL PAGE 19 
OF POOR QUALITY 

MEAN AXIAL VELOCITY PROFILES (CONT) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

n- 90° 
n = 270 » 

u = Q° 

« ~ 1 80 0 


FIG. 13 (CONT ' 


SYMBOL 

o 

• 

□ 

■ 

A 

▲ 

0 

♦ - 

RUN NOS 

14. 10, 42. 44 

15. 12. 

43. 45 

25. 26 

40 


i = 1 5? mm 


RUNS 14. 1-5. 25 i =-203 mm 


RUNS 10. 12. 2<> -l.i 41 




0 2 

0 4 

0 6 

0 3 ID 0O .02 

RADIUS RATIO, r/R 0 

0 4 

0 6 . 

0.8 1 

305 mm 



RUNS 42. 43 i - 400 mm 



RUNS 44. 45 




M 10 0 0 0 2 

RADIUS RATIO, r/R„ 


0 e 0.8 10 


83 - 3 - 76—29 


o: 













R83-915540-26 


ORIGINAL PAGE 59 
OF POOR QUALITY 


MEAN RADIAL VELOCITY PROFILES 

HORIZONTAL 1 VERTICAL^ 
TRAVERSE TRAVERSE 

OPEN SYMBOLS: 9 = 90* 9 a 0* 

SOLID SYMBOLS: 9 = 270* 9 = 180* 


SYMBOL 


RUN NOS. 


2 = 5 mm 


11,7,8.9 16,28,32,34 47,46,54 


RUNS 11. 16 


2 = 25 mm 


RUNS 7, 28. 47 



RADIUS RATIO, r/R 0 


z - 51 mm 

RUNS 8. 32, 46, 50. 51 z = 102 mm 

RUNS 9. 34.54 


RADIUS RATIO. r/R r 


83 - 3 - 76-28 


o 

















RADIAL VELOCITY, V — m/s 


o 


R83-91 5540-26 


ORIGINAL PAGE 19 
OF POOR QUALITY 

MEAN RADIAL VELOCITY PROFILES (CONT) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

# = 90* 

9 = 270° 

9 = 0° 

9 = 180? 


FIG. 14 (CONT.) 


z = 1 52 mm 


SYMBOL 


o • 0 


A A 


RUN NOS. 15,10.42,45 36.38 

RUNS 15.36 z = 203mm 


RUNS 10. 38. 41 


. * : » . . 

■ a S o'rCs'.A.' 11 


• A t g 

loo* 




® ® 0-2 M 0.6 - 0.8 1.0 O.O 0 2 0.4 0.6 0.8 10 

RADIUS RATIO, r/R„ 


z = 305 mm 


RUN 42 z =406 mm 


OS 1.0-00 

RADIUS RATIO, r/R r 


RUN 45 



2 04 . 06 0 


83 - 3 - 70-27 




o 








lGj 


R83-91 5540-26 


original page is 

OF POOR QUALITY 

MEAN AZIMUTHAL VELOCITY PROFILES 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

« = 90 ° 
n = 270° 

// = 0° 


SYMBOL 

0 • 

□ ■ 

A ▲ 

o ♦ 

RUN NOS. 

17, 3, 4, 5 

27, 31, 6 

47, 46, 33 

50, 54 

= 5 mm 


RUN 17 

z = 25 mm 



•o^ 

° V 


■'o* 

«D O • 


▼ 

51 

RUNS 3.27.47 



.8 8.2 8.4 8 6 8.e I 8 6.0 0.2 0.4 

RADIUS RATIO,. rJR n 


e.8 i. 


z = 51 mm 


RUNS 4, 31,46, 50. 51 


z = 1 02 mm 


RUNS 5. 6. 33. 54 


►|i<j***r*4 


18 1.1 II 1.2 

• 4 

It 

• 1 1.9 

RADIUS RATIO, r/R 0 



83 - 3 - 76-23 


0 







c 


o 


R83-91 5540i26 


ORIGINAL PAGE 19 
OF POOR QUALITY 

MEAN AZIMUTHAL VELOCITY PROFILES (CONT) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

fi = 90° 

0 = 270° 

0 = 0° 
0= 1Q0° 


FIG. 1 5 (CONT) 


z- 152 mm 


SYMBOL 


RUN NOS. 14,12.43,44 


A A 


RUNS 14,35 z — 203 mm 


RUNS 12. 3/. 40 


OCfc^ 





RADIUS RATIO, ,r/R 0 


z - 305 mm 

i — 

RUN 43 z = 406 mm 

""i r— — 

RUN 44 


OS 10 6 0 0 2 0 4 Of. . 0 . 

RADIUS RATIO, -i/Ba 


8 1 0 


83 - 3 - 76-26 


o 







RMS AXIAL VELOCITY FLUCTUATION, u' — m/s 


ORIGINAL PAGE fH 

R83-91 5540-26 OF POOR QUALITY FlG.1 

FLUCTUATING AXIAL VELOCITY PROFILES 



HORIZONTAL 

VERTICAL 


TRAVERSE 

traverse 

OPEN SYMBOLS: 

« = 90 ° 

Hs 0° 

SOLID SYMBOLS: 

H = 270° 

0 = 1 80 g 



RADIUS RATIO, r/R 0 


2 = 51 mm RUNS 4. 8. 23 z = 102 mm RUNS 5. 6. 9. 24 



&3-=3— 76— 25 






VELOCITY FLUCTUATION, u' — m/s 



















RMS RADIAL VELOCITY FLUCTUATION, v' — m/s 


R83-91 5540-26 


ORIGINAL PAGE 19 
OF POOR QUALITY 

FLUCTUATING RADIAL VELOCITY PROFILES 


FIG. 17. 



2 = 51 mm RUNS 8, 32, 46, 50, 51 2 = 102 mm RUNS 9. 34. 54 



0 0 0 2 0.4 0 6 e .8 1.0 0.0 0.2 1.4 16 0.8 1.8 

RADIUS RATIO, JiB 0 


83 - 3 - 76-22 









































AZIMUTHAL VELOCITY FLUCTUATION 





















INNER JET FLUID CONCENTRATION ALONG TEST SECTION CENTERLINE 



DISTANCE FROM TEST SECTION INLET, 






82 11 - 36-3 











mON FRACTION, f 


R83-91554Q-26 


FIG. 20 (CONT) 


ORIGINAL PAGE ffi 
OF POOR QUALITY 


MEAN INNER JET FLUID CONCENTRATION PROFILES (CONT) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

<9 = 90° 

0 = 270° 

0 = 0° 
6 * 180° 


SYMBOL O • □ 

RUN NOS. 25, 26 ; 


A | A 

36. 38 


z = 152 mm 


RUN 25, 35,36 z = 203 mm 



RUN 26. 37. 38 



RADIUS RATIO, r/R, 


82 - 11 - 36-4 









RMS INNER STREAM CONCENTRATION FLUCTUATION, 




R83-9 15540-26. 


ORIGINAL PAGE IS 
OF POOR QUALITY 

FLUCTUATING INNER JET FLUID CONCENTRATION PROFILES 


2 = 13 mrri 


r 

■* 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

0 = 90° 

0 = 270° 

0 = 0° 
0= 180° 


SYMBOL O 


A I A 0 ♦ 


RUN NOS. 20, 22, 23. 24 21 , 27. 31 , 33 29, 28, 32, 34 

RUNS 20. 21. 29, 30 ...z = 25mm 


RUNS 22. 27. 28 



0 0 0 2 0 4 


z = 51 mm 


-• • o ^ 

. A A 4 * ■ 

A ** A i 


a Q 








-1 NO 


R83-91 5540-26 


FIG. 21 (CONT) 


ORIGINAL PAGE 19 
OF POOR QUALITY 


FLUCTUATING INNER JET FLUID CONCENTRATION PROFILES (CONT) 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE 

OPEN SYMBOLS. 

(/ — 90 0 

0 = 0” 

SOLID SYMBOLS: 

0 = 270° - 

o 

O 

CO 


SYMBOL 

0 

• 

□ 

■ 

A 

▲ 

RUN NOS. 

25. 26 

35. 37 

36. 38 


z = 152 mm 


RUN 25. 35. 36 Z = 203 mm 


RUN 26. 37, 38 






o 


R83-915540-26 


2 = 5 rhm 


ORIGINAL PAGE IS 
OF POOR QUALITY „ 

MOMENTUM TRANSPORT RATE, uv, PROFILES 


HORIZONTAL 

VERTICAL 

TRAVERSE 

traverse 

6 = 90° 

11 = 0° 

0 = 270° 

ff = 180° 


SEE TABLE I FOR TRAVERSE DIRECTION 


SYMBOL 
RUN NOS. 


j 11,7.8.9 
RUNS 11. 16 



z =■ 25 mm 




•°o* 


0.2 

0.4 

0.6 

0 8 1.8 0.0 0 




RADIUS RATIO, r/R 0 

2 = 51 mm 



RUN 8 z = 102 mm 


U0 0 6 0 2 0 4 0 6 0 9 _ 1 fl 

RADIUS RATIO, r/R n 


82 10 33 13 


o 









MOMEN TUM TRANSPORT RATE, uv 


o 


R83-91 5540- 26 


ORIGINAL PAGE m 
OF POOR QUALITY 

MOMENTUM TRANSPORT RATE, uV, PROFILES (CONT.) 


z = 152 mm 


horizontal 

TRAVERSE 

VERTICAL 

TRAVERSE 

0 = 90° 

0 = 270 8 

0 = 0° 

0 = 180° 


o • 

□ i 

5. 10, 42, 45 

41 

RUN 15 2 

= 203 mm 


FIG. 22 (CONT.) 


RUNS 10. 41 




2 04 00 00 10 0.0 0.2. 0 4 0.6 0.3 1 

RADIUS RATIO, r/R n 


z = 305 mm 


RUN 42 z = 406 mm 


RUN 45 


RADIUS RATIO.. f/Rr 


83 — 3 — 76—18 










MOMENTUM TRANSPORT CORRELATION COEFFICIENT R 



82 — 10 - 33—19 






MOMENTUM TRANSPORT CORRELATION COEFFICIENT, 


R83-91 5540-26 


MOMENTUM TRANSPORT CORRELATION COEW&lflWJ-fiX, PROFILES (CONT. 


HORIZONTAL 

traverse 

VERTICAL 

traverse 

* = 90 e 
n = 270* 

8 = 0° 

8- 180° 


SYMBOL O • □ 

RUN NOS. 15, 10, 42, 45 


z = 152 mm 


RUN 15 — z = 203 mm 


3 0 00 



RUNS 10, 41 









MENTUM 


R63-91 5540-26 


ORIGINAL PAGE US 
OF POOR QUALITY 


z = 5 mrrl 


MOMENTUM TRANSPORT RATE, uw, PROFILES 


SYMBOL 
RUN NOS. 


HORIZONTAL 

VERTICAL 

TRAVERSE 

TRAVERSE 

0 = 90° 

o 

O 

: - 270° 

0 = 180° 


17,3,4,5 | 6 

RUN 17 z = 25 mm . 


68 e 2 8.4 8.6 0.8 1 8 08 02 04 8.6 88 J 

RADIUS RATIO, r/R 0 


z = 51 mm 


RUN 4 z = 102 mm 


RUNS 5. 6 


e i.e 0.0 o 2 

8 4 

8.6 

0 3 l_0 

RADIUS RATIO, r/R 0 



82 - 10 - 33-15 






MOMENTUM TRANSPORT RA 


ORIGINAL PAGE 19 
OF POOR QUALITY 


O 


R83-91 5540-26 


MOMENTUM TRANSPORT RATE, uw, PROFILES (CONT.) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

0 = 90' 
fl = 270' 

0 = 0' 
0=180' 


SYMBOL O • □ 

RUN NOS. 14, 12, 43, 44 


FIG. 24 (CONT' 


z - 1 52 mm 


RUN 14 z = 203 mm 


RUNS 12, 40 



o: 








MOMENTUM TRANSPORT CORRELATION COEFFICIENT 


o 


R83-91 5540-26. ORIGINAL PAGE S* 

OF POOR QUALITY 

MOMENTUM TRANSPORT CORRELATION COEFFICIENT, R uw , PROFILES 


HORIZONTAL 

VERTICAL 

TRAVERSE 

TRAVERSE 

OPEN SYMBOLS 

«s0° 

SOLID SYMBOLS ft -- e?u» 

V = 1 80 0 


2 -- 5 mm 


SYMBOL 
RUN NOS 


} /. 3. A , b 


RUN 17 2 = 25 Trim 


RADIUS RATIO, r/R 0 


z-51 mm 


2 — 1 02 mm 


RUNS 5. 6 


' 10 0 0 0.2: — 0 4 H • -0? 10 

RADIUS RATIO, r/R n 


82 10 33 21 


o 








MOMENTUM TRANSPORT CORRELATION COEFFICIENT 


ORIGINAL PAGE fc' 

OF POOR QUALITY F| _ ,,, rnNIT v 

B83-915540-26 FIG. 25 (CONT.) 

MOMENTUM TRANSPORT CORRELATION COEFFICIENT, R uw , PROFILES (CONT.) 


HORIZONTAL 

VERTICAL 

TRAVERSE 

TRAVERSE 

0 = 90° 

0 = 0° 

19 = 270° 

« 

o 

CO 

If 


SYMBOL O • □ ■ 

RUN NOS 14,12.43,44 40 


2 = 152 mm 


RUN 14 .... 2 = 203 mm 


RUNS 12. 40 



0.8 i .0 0 0 0 ,2 0.4 0.6 0.8 1.0 

RADIUS RATIO, r/R 0 


83 — 3 — 76—14 


o 








MOMENTUM TRANSPORT RATE, wv 


R83-9-1-5540- 2.6 


ORIGINAL PAGE (S 
OF POOR QUALITY 


MOMENTUM TRANSPORT RATE, wv, PROFILES 

I HORIZONTAL I VERTICAL 
TP AVERSE TRAVERSE 

OPEN SYMBOLS: «= 90° «~0" 

SOUP SYMBOLS « - 270° « •- 160” 


SYMBOL 


RUN NOS I 47. <16, 64 I 48* SO, 55* ~ 49* 51. 56* 



/ = 25 mm 


HUNS 4/ .18 49 


* DATA OBTAINED WITH 
OPTICAL AXIS A V DEG 
FROM TEST SECTION 
CENTERLINE 



e.e 02 e .4 0 6 0 8 i 


RADIUS RATIO,. ryR 0 


l - 51 mm 


RUNS 46. 50. 51, 52. 53 


7 = 102 mm 


RUNS 54. 55. 56 



RADIUS RATIO, .r/R. 


162 


H.1 1 58 1 















MOMENTUM TRANSPORT CORRELATION COEFFICIENT, R, 


o 


R83-915540-26 


ORIGINAL PAGE ES 
OF POOR QUALITY 

MOMENTUM TRANSPORT CORRELATION COEFFICIENT, R wv , PROFILES 



HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 



OPEN SYMBOLS: 

90° 

(1 = 0° 



SOLID SYMBOLS: 

a = 270° 

(1 = 1 80 6 



SYMBOL 


RUN NOS. I 47, 46. 54 


* DATA OBTAINED WITH 
OPTICAL AXIS ± 1 DEG 
FROM TEST SECTION 
CENTERLINE 


48* 50,55* 49* 51,56* 


z = 25 mm 



RUNS 47. 48. 49 



0.0 0 2 0.4 0 6 84 t 

RADIUS RATIO, r/R n 


z = 51 mm 


RUN 46, 50. 51, 52. 53 z = 102mm 


RUNS 54. 55. 56 




RADIUS RA T !0, r/R, 


83—1 — 68—2 


G 















HORIZONTAL 

TRAVERSE 

OPEN SYMBOLS: « = 90 c 

SOLID SYMBOLS: « = 270° 


VERTICAL 

TRAVERSE 

« = 0 ° 

«= 180° 

















HORIZONTAL VERTICAL 
TRAVERSE TRAVERSE 


OPEN SYMBOLS: tf = 90° 

SOLID SYMBOLS: « = 270" 


il. - 0° 
li= 180° 











RADIAL MASS TRANSPORT CORRELATION 
COEFFICIENT, R wf 


R83-91 5540-26 


FIG. 29 (CONT.) 


ORIGINAL PAGE 
OF POOR QUALITY 


RADIAL MASS TRANSPORT CORRELATION COEFFICIENT, R vf , PROFILES (CONT.) 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE 

OPEN SYMBOLS: 

u = 90° 

u- 0° 

SOLID SYMBOLS: 

II = 270° 

it 

CD 

O 

o 


z =.152 mm 

0 4 


hUN 36 z = 203 mm 


Dl IM od 





















R83-91 5540-26 


ORIGINAL PAGE 19 
OF POOR QUALITY 


EIG. 30 (CONI.) 


AXIAL MASS TRANSPORT RATE, u?, PROFILES (CONT.) 



HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

OPEN SYMBOLS. 
SOLID SYMBOLS: 

0 = 90° 

61 = 270° 

0 = 180' 


0.04 


z = 51 mm 


RUN 23 


i - 102 mm 


RUN 24 








TRANSPORT CORRELATION COEFFICIENT, R 
















AXIAL MASS TRANSPORT CORRELATION 
COEFFICIENT, R llf 


R83?91554Q-26 


FIG, 31 (CONT.) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


AXIAL MASS TRANSPORT CORRELATION COEFFICIENT, R uf , PROFILES (CONT.) 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE 

OPEN SYMBOLS: 

a = 90° 

fl = 0° 

SOLID SYMBOLS: 

o 

O 

f"- 

OJ 

II 

o 

O 

CO 

1! 



RADIUS RATIO, r/R 0 


83—2-55- 1 tB 





AZIMUTHAL MASS TRANSPORT RATE, wf — m/s 


R83-91 5540-26 


ORIGINAL PAGE 19 
OF POOR QUALITY 


AZIMUTHAL MASS TRANSPORT RATE, wf, PROFILES 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

= 90° 

0 = 270° 

0 = 0 ° 

0 = 180° 


RUN 29 z = 25 mm 



RUN 27 



0.6 0.S 1.OO.0 0.2 0 4 

RADIUS RATIO, t/B 0 


0 6 OS 


RUN 31 z = 102 mm 


RUN 33 



RADIUS RATIO, r/R r 


83-2-55— 23A 


o; 








AZIMUTHAL MASS TRANSPORT RATE, 


R83-91 5540-26 


FIG. 32 (CONT.) 


ORIGINAL PAGE 19 
OF POOR QUALITY 


AZIMUTHAL MASS TRANSPORT RATE, wf, PROFILES (CONT.) 



83-2-55-23B 





AZIMUTHAL MASS TRANSPORT CORRELATION! COEFFICIENT 












AZIMUTHAL MASS TRANSPORT 
CORRELATION COEFFICIENT, R,. 


R83-9 15540- 26 


FIG. 33 (CONT.) 


ORIGINAL PAGE 19 - 
OF POOR QUALITY 


AZIMUTHAL MASS TRANSPORT CORRELATION COEFFICIENT, R w j, PROFILES (CONT.) 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE 

OPEN SYMBOLS. 

« = 90° 

« 0° 

SOLID SYMBOLS: 

"-■ 270 ° 

«= 180° 



RADIOS. RATIO. r/R Q 


83—2 -55—1 3B 














R83-91 5540-26 


ORIGINAL PAGE IS 
OF POOR QUALITY 

RADIAL VELOCITY PROBABILITY DENSITY FUNCTIONS 

AXIAL LOCATION: 25 mm — 


AVERAGES OF DATA FROM RUNS 7 (O) AND 28 (□) 



,R o = 


Av = 0.1 m/s 


































R83-91 5540-26 


ORIGINAL PA.GE IS 
OF POOR QUALITY 

SKEWNESS OF AXIAL VELOCITIES PROFILES 


FIG. 3 


SYMBOL 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

« s 90° 

« = 270“ 

= 0° 

0=180°- 


O • □ 


^ I A 0 ♦ 


RUN NOS 11,3,4,5 16,7,8,6 17,22,23,9 


z = 5 mm 


RUNS 11, 16, 17 2 = 25 min 


RUNS 3. 7. 22 



RADIUS RATIO, r/R Q 


2 = 51 mm 

i 

RUNS 4, 8. 23 — 2 = 102 mm 

| 1 " 

RUNS 5. 6. 9. 24 




RADIUS RATIO. r/R 0 

— 


83 - 3 - 76-11 









o 


R83-91 5540- 26 


ORIGINAL PAGE fa 
OF POOR QUALITY 

SKEWNESS OF AXIAL VELOCITIES PROFILES (CONT.) 


FIG 38 (CONT.) 


horizontal 

TRAVERSE 

VERTICAL 

TRAVERSE 

// = 90 0 
0 == 270° 

// s 0° 

1 80° 


Symbol I o 


^ a o nr 


RUN N0 $ 14.10.42.44 15.-12, 43. 45 25~26 | 

4~152mm RUNS 14, 15. 25 z = 203 mm 


RUNS 10. 12. 26. 40 41 



RADIUS RATIO. r/R Q 


z = 305 mm 



RUNS 42, 43 2 = 406 mm 

— r | 

RUNS 44. 45 


>8 1 0 0.0 0.2 0.4 

0.6 0.0 1.0 

RADIUS RATIO, r/R 0 

83-3- 76-11 A 















RS3-9 15540* 26 



ORIGINAL PAGE 19 
OF POOR QUALITY 

KURTOSIS OF AXIAL VELOCITIES PROFILES (CONT.) 


OPEN SYMBOLS 


horizontal^ 

traverse 

VERTICAL 

traverse 

// = 90°~ 
// = 270° 

H = 0 6 
«=180° 


FIG. 39 (CONT.) 


_ SY »». | o I . d[ T | A I 4 I O I T 

WNNOS | 14. 10. 42, 44 15 . ,2. ,3 „ ^ ± 


'* - 1 52 mm 


RUNS 14. 15.25 2 = 203 mm 


RUNS 10. 12. 26. 40. 41 


-a 4| .4 A 


4 ,.*» 


o • 


«* I !, .*8 V ; *■:„* 



RADIUS RATIO, r/R. 


83-3-76 -12A 


o: 






083 * 91 5540 - 26 


Tpoob quality 

SKEWNESS OF RADIAL VELOCITIES PROFILES 




SYMBOL 


HORIZONTAL VERTICAL 
TRAVERSE TRAVERSE 

OPEN. SYMBOLS « = 90° //:.£)» 

SOLID SYMBOLS </ = 270° '/=180° 


o • A A o ♦ V * 


RUN NOS 11.7.8.9 16.28,32,34 47,46,54 


7 - 5 mm 


RUNS 11. 16 7 = 25 mm 


RUNS 7. 28. 47 



RADIUS RATIO, r/R Q 


z - 51 mm 

RUNS 8. 32. 46. 50, 51 z -102 mm 

RUNS 9. 34. 54 


RADIUS RATIO, r/R r 


83 - 3 - 76-10 


o 








R83*915540=26 


ORIGINAL PAGE ISI 
OF POOR QUALITY 


SKEWNESS OF RADIAL VELOCITIES PROFILES (CONT.) 



HORIZONTAL 

VERTICAL 


traverse 

traverse 

OPEN SYMBOLS 

(/~90 c 

ft = 0° 

SOLID SYMBOLS 

ft - 270 6 

180° 


FIG. 40 (CONT.) 


z = 1 52 mm 


Symbol o • qJ J T 

RUN N0S 15.10.42 36,38.45 

__ RUNS 15,36 2 = 203 mm 


A A 


RUNS 10. 38. 41 



RADIUS RATIO, r/R Q 

z = 305 mm 

run 42 2 = 406 mm 


1 0 0.0 0.2 

RADIUS RATIO, r/R 0 

0 4 

0 6 

0.e i.0 




03 - 3- 76— 1 0 A 


o: 



















KURTOSIS 


c 


m 


R83-91 5540-26 


ORIGINAL PAGE 19 
OF POOR QUALITY 

KURTOSIS OF RADIAL VELOCITIES PROFILES (CONT.) 


FIG. 41 (CONT.) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

H ~ 90° 
H = 210° 

u = 0° 

H 1-180' 


z = 152 mm 


SYMBOL 
RUN NOS 


15,10,42 
RUNS 15, 36 


□ ■ 

36. 38, 45 


z - 203 mm 


RUNS 10, 38. 41 


I Q 8 


□ 00 ■■ 
0 D D 


0, B/ 0jt 



RADIUS RATIO, r/R Q 


z = 305 mm 

RUN 42 z = 406 mm 

RUN 45 


RADIUS RATIO, r/R> 


83— 3—76— 9A 


o: 































ORIGINAL PAGE 19 
OF POOR QUALITY 


o 


R83-91 5540-26 ur ruwn FlG .42 (CONT’ 

SKEWNESS OF AZIMUTHAL VELOCITIES PROFILES (CONT.) 



HORIZONTAL 

VERTICAL 


TRAVERSE 

traverse 

OPEN SYMBOLS 

« = 90“ 


SOLID SYMBOLS 

« = 27 0° 

o — 160 * 


SYMBOL 


0 

• 

□ 

■ 

A 

A 

14, 12, 

43. 44 

35. 

37 

40 


z = 152 mm 


RUNS 14,35 z - 203 mm 



RUNS 12. 37. 40 



.0 0.2 0 

.4 0.6 

0.8 1.0 0 . 

0 .. 0.2 0 

.4 0.6 

0.8 1 . 



RADIUS RATIO , r / R Q 



z = 305 mm 

r —-- ■ 1 


RUN 43 

z = 406 mm 


RUN 44 




RADIUS RATIO, r/R, 


83— 3—76-8A 


O 








URTOSIS, K. 


c: 


o 


R83-91 554CF26 


ORIGINAL. PAC£ K 
OF POOR QUALITY 

KURTOSIS OF AZIMUTHAL VELOCITIES PROFILES 


HORIZONTAL 

traverse 

VERTICAL 

TRAVERSE 

«-j=90° 
ft = '270 - 

«=0° 

'/ = 180° 


FIG. 43 


SYMBOL 
RUN NOS 
/. - 5 mm 


IV, 3. 4. 5 


27. 31, 6 
RUN 17 


Li_ 

▲ 

0 

♦ | 

47. 46. 33 

50. 54 


7 - 2h mm 


RUNS 3. 27. 4 7 


0 *>• 






* D 0 


• ¥ oo#° ^ ;•* • • 

A * ■ 


• » 8 » 


e 8 0.2 0 4 0 6 « 8 1 0 0.0 0.2 0.4 0.6 0.8 

RADIUS RATIO. r/R 0 


2 = 51 mm 


RUNS 4,-31. 46. 50. 51 


2 = 102 mm 



RUNS 5. 6. 33. 54 


RADIUS RATIO. r/R, 


83 3 76 7 


o: 








R83-91 5540-26 


ORIGINAL PAGE 12J 
0F POOP QUALITY 

KURTOSIS OF AZIMUTHAL VELOCITIES PROFILES (CONT.) 


FIG. 43 (CONT.) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

ft » 90° 
a = 270* ' 

« 0® 
(is 180® 


z = 1 52 mm 


SYMBOL 


O • □ ■ A A 


RUN NOS 14,12,43,44 35,37 

RUNS 14, 35 z = 203 mm 


RUNS 12. 37. 40 


t * / * H 





RADIUS RATIO, r/R 0 


z = 305 mm 

RUN 43 z = 406 mm 

RUN 44 


® fi o # 0 o 

• <>•© ° O 


I . I . I . I . I . I . I . I . I . 


RADIUS RATIO, r/R, 


83— 3—76— 7A 












R83-9 15540- 26 


ORIGINAL PAGE IU 
OF POOR QUALITY 


SKEWNESS OF INNER JET FLUID CONCENTRATION PROFILES 


SYMBOL 


2 -- 1 3 mm 


HORIZONTAL 

traverse 

VERTICAL 

traverse 

0= 90° 

0 = 270° 

« = 0° 
II- 180° 


RUN NOS | 20. 22. 23. 24 21 , 27. 31 . 33 I 29. 28. 32. 3- 


RUNS 20. 21. 29. 30 


z = 25 mm 


RUNS 22. 27. 28 


2 = 51 mm 


RADIUS RATIO. r/R 

RUNS 23. 31. 32 2 = 1 


'"‘O 

2 = 102 mm 


RUNS 24. 33. 34 


** ** * 


-4 

V 2 

a 

D 

.at D I ? If *>fi* 


RADIUS RATIO, r/R r 


82 11 - 36-7 


















R83-91 5540-26 


FIG. 44 (CONI.) 


ORIGINAL PACE t® 
OF POOR QUALITY 


SKEWNESS OF INNER JET FLUID CONCENTRATION PROFILES (CONT) 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE 

OPEN SYMBOLS: 

a=90° 

0 = 0° 

SOLID SYMBOLS: 

0 = 270° 

0*180°. 


SYMBOL 

o 


° 

■ 

A 

A 

RUN NOS. 

25. 26 

35. 37 

— 

36. 

3ft 


15 


z = 152 nnm 


RUN 25. 35. 36 


z = 203 mm 


RUN 26. 37. 38 



























o 





© 








R83-915540-26 


FIG. 46 


ORIGINAL PAGE Si 
OF POOR QUALITY 


TURBULENT MOMENT UM TRANS PORT RATE, u*, PROBABILITY DENSITY FUNCTIONS 

AXIAL LOCATION: 25 mm 


DATA TROM RUN 7 


r/R Q 

uv 

°uv 

> 

D 

CO 

c 

< 

0.00 
0.13 
0.24 
0.36 
0.62 ... 
0.68 

-0.0003 

0.0118 

0.0083 

*0.0623 

*0.0119 

0.0254 

0.060 
0.055 
0.072 
0.164 
0 041 
0.085 

1.58 

*0.04 

0.25 

-1.11 

-1.44 

1.26 

19.27 
6.19 
10.15 
8.21 
10.42 
12.50 - 


A(uv) = 0.02 m 2 /s 2 




[/R o = 0.24. 


df 


°o 


ax85$? 1 fooaQnl -a 


2/ e 2 


MOMENTUM TRANSPORT RATE, uv — mWs 




0.5 -0.4 -0.3 -0.2 -0.1 0. 0.1 Q.2 0.3 -0.3 -0. 2:0.1 -0 0.1 0.2 0.3 



MOMENTUM TRANSPORT RATE, uv — m 2 /s 2 


83 719 7 






























o 


R83-91 5540-26 


OF POOR QUALITY 


FIG. 4 


TURBULENT RADIAL MASS TRANSPORT RATE, vf, PROBABILITY DENSITY FUNCTIONS 

AXIAL LOCATION: 25 mm 


data From run 28 


•0.0001 0.018 

0.0002 0.018 

0.0004 0.056 

0.0041 0.068 

•0.0004 0.003 

•0.0001 0.002 

A(vf) = 0.04 m/S 


f/R o = 0.0 


r/R 0 = 0.13 


r/R o = 0.24 





RADIAL MASS TRANSPORT RATE, vf — m/s 





•0.3 -0.2 -0.1 0 0.1 0.2 0.3 


-0.3 -0.2 -0.1 0 0.1 0.2 0.3 


-0.a.-0.2 -C.1 0 0.1 0.2 0.3 


RADIAL MASS TRANSPORT RATE, vf — m/s 


83 - 7 - 19-4 





R83-91 5540-26 


ORIGINAL PAG£ ft 
OF POOR QUALITY 


FIG. 


TURBULENT AXIAL MASS TRANSPORT RATE, uf, PROBABILITY DENSITY FUNCTIONS 

AXIAL LOCATION: 25-rnm 


Data from run 22 


r/R 0 

uf 

°ul 

S ul 

K u1 

0,00 

0.13 

0.24 

0.36 

0.52 

0.68 

0.0052 

0.0137 

0.0318 

*0.0020 

0.0001- 

*0.0005 

0.039 
0.089 . 
' 0.10S 

0.045 
0.003 
0.003 

9.17 
1.13 
0.82 ... 
-0.38 
6.30 
*1.65 

110.6 

21.7 

6.9 

6.6 

137.7 

18.4 


*i(uf) = 0.04 m/s 



AXIAL MASS TRANSPORT RATE, uf — m/s 


83 - 7 — 19 — 3 







R83-91 5540-26 


ORIGINAL. PAGE fS 
OF POOR QUALITY 


TURBULENT AZIMUTHAL MASS TRANSPORT RATE, wl, 
PROBABILITY DENSITY FUNCTIONS 

AXIAL LOCATION; 25 mm 


DATA FROM RUN 27 


r/R 0 

wf 

a wf 

s wl 

K wl 

0.0 

-0.0001 

0.022 

•1.56 

42.2 

0.13 

0.0002 

0.028 

•1.30 

18.2 

0.24 

0.0187 

0.057 

1.39 

15.0 

0.36 

*0.0005 

0.041 

0.76 

12.8 

0.52 

0.0003 

0.006 

3.76 

75.1 

0.68 - 

0.0003 

0.003 

•2.78 

30.2 


A(wf) = 0.04 m/s 


AZIMUTHAL 


o 




RAL MOMENT, 


R03-91 5540-26 


ORIGINAL PAGE IS 
OF POOR QUALITY 


SECOND CENTRAL MOMENT OF uv TURBULENT TRANSPORT RATE PROFILES 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE' 

OPEN SYMBOLS: 

0 = 90° 

0 = 0° 

SOLID SYMBOLS' 

270° 

ii 

CD 

o 

© 



RADIUS RATIO, r/R 0 








SECOND CENTRAL MOMENT 










R83-91 5540-26 


ORIGINAL PAGE IS 
OF POOR QUALITY 


SKEWNESS OF uv TURBULENT MOMENTUM TRANSPORT RATE PROFILES 



HORIZONTAL 

VERTICAL 


TRAV6RSE 

TRAVERSE 

OPEN SYMBOLS 

n a 90 0 

i/ = 0° 

SOLID SYMBOuS 

ii - 270‘ 

»= 180 c 


SYMBOL 

0 

• 

□ 

■J 

RUN NOS 

11.7. 8. 9 

1(3 


FIG. ! 


? - b mm 


RUNS 11, 16 


/ -- 95 mm 


RUN / 







ORIGINAL PAGE ES 
OF POOR QUALITY 


o 


R83*91 5540*26 FJG. 53 (CONT' 

SKEWNESS OF Qv TURBULENT MOMENTUM TRANSPORT RATE PROFILES (CONT.) 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

0 = 270° 

0 s 180° 


SYMBOL O • 

RUN NOS 15.10.42,45 


4-15? mm 


RUN 15 z — 203 mm 


RUNS 10. 41 



RADIUS RATIO, r/R Q 


z - 305 mm 

RUN 42 2 = 406 mm 

RUN 45 


RADIUS RATIO, r/R, 


83-3 76-6A 


o 








ORIGINAL PAGE IS 
OF POOR QUALITY 


R03-91 5540* 26 OF POO * ■ ■ 

KURTOSIS OF uv TURBULENT MOMENTUM TRANSPORT RATE PROFILES 


SYMBOL 


HORIZONTAL 

TRAVERSE 

VERTICAL 

TRAVERSE 

H a 90 0 

(t = 270° 

= o 0 

ff » 1 80 c 


o • 

□ i 

11. h 8. 9 

1G 

1.16 z 

mm 



RADIUS RATIO, r/R 0 

/. - 51 mm 

RUN 8 / = 102 mm 


RADIUS RATIO, r/R r 


83-3 76 5 



















R83-91 5540-26 


ORIGINAL PAGE 19 
OF POOR QUALITY 


F 


SECOND CENTRAL MOMENT OF uw TURBULENT TRANSPORT RATE PROFILES 


HORIZONTAL 

VERTICAL 

TRAVERSE 

traverse 

0 = 90° 

0 = 0° 

0 = 270° 

0=180° 


SYMBOL 
RUN NOS. 


2 = 5 mm 


17 , 3 , 4 , 5 


RUN 17 Z = 2b mm 



10 0 8 0 2 

RADIUS RATIO, r/R 0 


0 8 i . 0 


82 10 33 - 1 7 


O 







SECOND CENTRAL MOMENT, 












o 


ORIGINAL PAGE W 

R83-91 5540-26 OF POOR QUALITY 

SKEWNESS OF uw TURBULENT MOMENTUM TRANSPORT RATE PROFILES 


HORIZONTAL 

traverse 

VERTICAL - 
TRAVERSE 

a = 90 ° 

0*2 70" 

u = 0‘ 
U= 180' 


SYMBOL 
RUN NOS 


O • 

17. 3, 4, 5 


z = 5 mm 


RUN 17 2 = 25 mm 



RADIUS RATIO, r/R 0 


z - 51 mm 

f 

RUN 4 2 = 1 02 mm 

— — - ■ i f 

RUNS 5. 6 


RADIUS RATIO, r/R r 


83-3 76-4 


o; 









ORIGINAL PAGE IS 
OF POOR QUALITY 


o 


R83-91 5540-26 . FIG, 56 (CONT.) 

SKEWNESS.QFJR5 TURBULENT MOMENTUM TRANSPORT RATE PROFILES (CONT t ) _ 


HORIZONTAL 

TRAVERSE 

VERTICAL-! 
TRAVERSE ' 

ft = 90° 
ft = 2T0? 

0 = 0* 
0*180° 


SYMBOL 
RUN NOS 


14. 12. 43. 44 


z~ 152 mm 


7. - 203 mm 


RUNS 12. 40 



RADIUS RATIO, r/R 0 


z = 305 mm 

RUN 43 t = 406 mm 

RUN 44 


RADIUS RATIO, r/R f 


83— 3—76— 4A 


o 








KURTOSIS. K..-7- 


O 


ORIGINAL PAGE IS 

R83-91 5540-26 .... OF POOR QUALITY F4Gt g 7 

KURTOSIS OF uW TURBULENT MOMENTUM TRANSPORT RATE PROFILES 



horizontal 

VERTICAL 


traverse 

TRAVERSE 

OPEN SYMBOLS 

( 1 * 90 ° 

//SO 8 

Solid symbols. 

o 

o 

C\i 

II 

4X5 

o 

o 

CO 

u 


Symbol o • □ a 

RUN NOS 17.3,4.5 6 


2 z = 5 mm 


RUN 17 z = 25 mm 


• ° O 


o o 


°? Vo® *°. • 

o 

0 # O o 


• o o 


• • • 

°* O** 


t ^ • 


l 0 0 82 04 lU 88 10 00 0 2 .0 4 0 6 0 8 1 

RADIUS RATIO, x/R n 


z = 51 mm 



2-102 mm 


RUNS f) b 


• 0 » 

• • . 


• 2 ■ • O 


0 0 2 0 4 0 e 0 8 1-0 0 0 0 2 

RADIUS RATIO, r/R 0 


212 


o 





o 


QRIG5NAL FAGt V'o 

. Op POOR QUALITY 

R83-91 5540-26- , 

SECOND CENTRAL MOMENT OF wv TURBULENT TRANSPORT RATE PROFILES 



HORIZONTAL 

VERTICAL 


TRAVERSE 

Traverse 

OPEN SYMBOLS: 

90° 


SOLID SYMBOLS' 

I'D 

-nJ 

O 

o 

0 ~ \ 00 ° 


SYMBOL 


RUN NOS; I 47. 46. 54 


48,50,55 49,51,56 


1 = 25 mm 



RUNS 4/, 48. 49 


0 • A 


0 0 0.2 0 4 0 6 18 | 

RADIUS RATIO, r/R n 


z = 51 mm 


RUNS 4G, 50. 51. 52. 53 z = 102 i 


RUNS 54. 55. 56 


til |oH°* “••VfV* • 8o 3 ° 3 

i l -t - 1 i 1 i 1 » I i 1, i 1 i 1 i I i I i I r . i . i . i , i , i . i . i , i . 
0 0 0 2 0 4 If - If 1 0 0 0 0 2 0 4 0 6 0 8 

RADIUS RATIO, r/R Q 


PRECEDING PAGE BLANK NOT FTLMFn 


83 1 58-3 


O 



















original PAGE Ift 
OF POOR-QUALITY 


O 


R83-915540-26 

SKEWNESS OF wvTURBULENT MOMENTUM TRANSPORT RATE PROFILES 


HORIZONTAL 

TRAVERSE 

VERTICAL’ 

TRAVERSE 

H = 90* 

U = i 70° 

li = 180° 




SYMBOL O 


HUN NOS 47,46.54 48,50,55 I 49,51,56 


z = 25 mm 


RUNS 47. 48. 49 


ft 4 ° «* 

^ BAP A 


0 0 0.2 0.4 9.6 9.9 1 

RADIUS RATIO, r/R n 


2 -- 51 mm 


RUNS 46, 50, 51. 52, 53 z = 102mm 


RUNS 54, 55. 56 


0.9 1.9 0.9 0.2 

RADIUS RATIO, r/R d 

9 4 

9 6 

0.9 1.0 




83 - 3 - 76-2 


G 

















KURTOSIS, K, 


R83-915540-26 


ORIGINAL PAGE IS 
OF POOR QUALITY 


KURTOSIS OF. wv TURBULENT MOMENTUM TRANSPORT RATE PROFILES 


SYMBOL O • □ 

RUN NOS 47. 46. 64 48 


HORIZONTAL 

TRAVERSE 

vertical' 

TRAVERSE 

« = 90 # 
<t~27Q° 

it ~ o e 
« = 180° 


■ ^ 

A < 

5 49. St 

56 


z-25 mm 


RUNS 4 7. 48 46 


• o* 4 

■* * 0 

O a 

18 r o»o 0 v a □ 

4 • A 


00 0 2 0 4 0.6 0.8 i 

RADIUS RATIO, r/R rt 


z = 51 mm 


RUNS 46; 50.51.52.53 


2 = 1 02 mm 


RUNS 54. 55. 56 


■•« f! j i°« 

rt a * 


o o . ♦ . 

% $ 
*o o 


- 4 0 A A 

* fo 0 6 § ° 


* 8 8 2 0 4 0 0 • • t • • 0 0 0.2 0 4 0 6 

RADIUS RATIO. r/R n 


83 - 3 - 2 ( 5 - 1 


o 























R83-91 5540- 26 


ORIGINAL PA£»£ 

OF POOR QUALITY 

SKEWNESS OF TURBULENT RADIAL MASS TRANSPORT RATE PROFILES 



HORIZONTAL 

VERTICAL 


TRAVERSE 

TRAVERSE 

OPEN SYMBOLS: 

0-90° 

<l~Q° 

SOLID SYMBOLS: 

0 = 270° 

o 

O 

CO 

II 


f. - 1 3 mm 


RUN 30 


z = 25 mm 


RUN 28 








R83-915540-26 


ORIGINAL PAGE fS 
OF POOR QUALITY 


FIG. 62 (CONT. 


SKEWNESS OF TURBULENT RADIAL MASS TRANSPORT RATE PROFILES (CONT.) 


HORIZONTAL 

VERTICAL 

TRAVERSE 

TRAVERSE 

« = 90° 

» = 0° 

« = 270° 



z - 1 52 mm 


RUN 36 z *= 203 mm 


RUN 38 


1 10 0.0 0,2 0 4 0,6 0.9 i.e 

RADIUS RATIO, r/R n 


83-2 -55-lflU 


o 





KURTOSIS 












R83-91 55*0,26 


ORIGINAL PAGE fS 
OF POOR QUALITY. 


FIG. 63 (CONT.; 


KUBIOSIS OF TURBULENT RADIAL MASS TRANSPORT .RATE PROFILES (CONT.) 
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SECOND-CENTRAL .MOMENT OF TURBULENT AXIAL MASS TRANSPORT 
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SKEWNESS OF TURBULENT AXIAL MASS TRANSPORT RATE PROFILES 
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KURTOSIS OF TURBULENT AXIAL MASS TRANSPORT RATE PROFILES 
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KURTOSIS OF TURBULENT AZIMUTHAL MASS TRANSPORT RATE PROFILES (CONT.) 
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MEAN RADIAL VELOCITY PROFILES 

SWIRLER EXIT PLANE AT z = Omm OATA OBTAINED AT z = 5 mm SEE FIG. 74 FOR VANE ORIENTATION 
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FLUCTUATING RADIAL VELOCITY PROFILES 

SWIRLER EXIT PLANE AT z = 0 mm DATA OBTAINED AT z - 5 mm SEE FIG. 74 FOR VANE ORIENTATION 
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SWIRLER EXIT PLANE AT 2 -=-51 mm DATA OBTAINED AT 2 = 5 mm SEE FIG. 74 FOR VANE ORIENTATION- 
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